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Abstract A tutorial introduction is given to q-special functions and to q-analogues of the 
classical orthogonal polynomials, up to the level of Askey- Wilson polynomials. 



0. Introduction 

It is the purpose of this paper to give a tutorial introduction to g-hypergeometric functions 
and to orthogonal polynomials expressible in terms of such functions. An earlier version 
of this paper was written for an intensive course on special functions aimed at Dutch 
graduate students, and later it was part of the lecture notes of my course on "Quantum 
groups and (/-special functions" at the European School of Group Theory 1993, Trento, 
Italy. Eventually, this paper may appear as part of a book containing the notes of the 
main courses at this School. 

I now describe the various sections in some more detail. Section 1 gives an introduction 
to g-hypergeometric functions. The more elementary (/-special functions like (/-exponential 
and g-binomial series are treated in a rather self-contained way, but for the higher q- 
hypergeometric functions some identities are given without proof. The reader is referred, 
for instance, to the encyclopedic treatise by Gasper & Rahman [15]. Hopefully, this section 
succeeds to give the reader some feeling for the subject and some impression of general 
techniques and ideas. 

Section 2 gives an overview of the classical orthogonal polynomials, where "classical" 
now means "up to the level of Askey- Wilson polynomials" [8] . The section starts with the 
"very classical" situation of Jacobi, Laguerre and Hermite polynomials and next discusses 
the Askey tableau of classical orthogonal polynomials (still for q = 1). Then the example of 
big q- Jacobi polynomials is worked out in detail, as a demonstration how the main formulas 
in this area can be neatly derived. The section continues with the g-Hahn tableau and 
then gives a self-contained introduction to the Askey- Wilson polynomials. Both sections 
conclude with some exercises. 

Notation The notations N := {1, 2, . . .} and Z + := {0, 1,2,.. .} will be used. 

Acknowledgement I thank Rene Swarttouw for commenting on a preliminary version of 
this paper. 
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1. Basic hypergeometric functions 

A standard reference for this section is the recent book by Gasper and Rahman [15]. In 
particular, the present section is quite parallel to their introductory Chapter 1. See also 
the useful compendia of formulas in the Appendices to that book. The foreword to [15] by 
R. Askey gives a succinct historical introduction to the subject. For first reading on the 
subject I can also recommend Andrews [2] . 

1.1. Preliminaries. We start with briefly recalling the definition of the general hyper- 
geometric series (see Erdelyi e.a. [13, Ch. 4] or Bailey [11, Ch. 2]). For a G C the shifted 
factorial or Pochhammer symbol is defined by (a)o := 1 and 



(a) k :=a(a+l)...(a + k-l), k = 1,2, 



The general hypergeometric series is defined by 



r- 1 s 



ai, . . . , a r 



E Mk ■ ■ ■ Mk k / -j -j \ 

, =0 (6i)*...(6.)**l ' ' 



Here r, s G Z + and the upper parameters a\, . . . , a r , the lower parameters b\, . . . , b s and 
the argument z are in C. However, in order to avoid zeros in the denominator, we require 
that b\, . . . , b s ^ {0, —1, —2, . . .}. If, for some % = 1, . . . , r, is a non-positive integer then 
the series (1.1) is terminating. Otherwise, we have an infinite power series with radius of 
convergence equal to 0, 1 or oo according to whether r — s — 1>0, = 0or<0, respectively. 
On the right hand side of (1.1) we have that 

(k + l)th term {k + a±) ... (k + a r ) z 

kth term ~ (k + h) ... (k + b s ) (k + 1) ^ ' ' 

is rational in k. Conversely, any rational function in k can be written in the form of the 
right hand side of (1.2). Hence, any series J2T=o Ck with cq = 1 and Ck+i/ck rational in k 
is of the form of a hypergeometric series (1.1). 

The cases o-^o and iF are elementary: exponential resp. binomial series. The case 
2-F1 is the familiar Gaussian hypergeometric series, cf. [13, Ch. 2]. 

We next give the definition of basic hypergeometric series. For a, q G C define the 
q-shifted factorial by (a;q)o := 1 and 

(a; q) k := (1 - a)(l - aq) . . . (1 - ag*" 1 ), k = 1, 2, . . . . 

For \q\ < 1 put 

oo 

(a; q)oc := Y[(l-aq k ). 

fc=0 

We also write 



(a 1 ,a 2 ,...,a r ]q) k := (a 1 ;q) k (a 2 ;q)k ■■■ (a r ;q)k, k = 0,1,2,... or oo. 
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Then a basic hypergeometric series or q-hypergeometric series is defined by 
ai, . . . , a r 



-<f>s 



h,...,b. 



= r s (ai,...,a r ;6i,...,6 s ;q',2;) 

(ai,...,a r ; ? ) fc j ( _ 1)fcg * ( *_i )/2 ji+.-r r?sGZ+ . 



fc=0 



(6i,...,6 a ,g; 



(1.3) 



On the right hand side of (1.3) we have that 

(Jfe + l)th term (1 - ai q k ) ... (1 - a r q k ) (-q k ) 1 + s ~ r z 



(1.4) 



fcth term (1 - b x q k ) ... (1 - 6 s g fc ) (1 - q^ 1 ) 

is rational in q k . Conversely, any rational function in q k can be written in the form of the 
right hand side of (1.4). Hence, any series X^fcLo ° k with cq = 1 and Ck+i/ck rational in q k 
is of the form of a (/-hypergeometric series (1.3). This characterization is one explanation 
why we allow q raised to a power quadratic in k in (1.3). 
Because of the easily verified relation 

(a; = (-!)* o fc «-* (fc - 1)/2 (a" 1 ;?),, 

any series (1.3) can be transformed into a series with base q . Hence, it is sufficient to 
study series (1.3) with \q\ < 1. The tricky case \q\ = 1 has not yet been studied much and 
will be ignored by us. Therefore we may assume that \q\ < 1. In fact, for convenience we 
will always assume that < q < 1, unless it is otherwise stated. 
In order to have a well-defined series (1.3), we require that 

^ l,5 _1 ,g" 2 ,... . 

The series (1.3) will terminate iff, for some i = l,...,r, we have G {1, q~ 2 , . . .}. 
If ai = q~ n (n = 0, 1, 2, . . .) then all terms in the series with k > n will vanish. In the 
non-vanishing case, the convergence behaviour of (1.3) can be read off from (1.4) by use 
of the ratio test. We conclude: 



convergence radius of (1.3) 



00 if r < s + 1, 

1 if r = s + 1, 
if r > s + 1. 



We can view the (/-shifted factorial as a (/-analogue of the shifted factorial by the limit 
formula 

= (a) fc := o(o + 1) • • • (a + k - 1). 



lim 



Hence 



(1-g)* 

is a (/-analogue of r F s by the formal (termwise) limit 



lim, 

9ti 



7 ai 



(/ b S ...,(/" 



l+s— r , 





«1, • 




— p 

r- 1 s 






.61,. 


h ]Z 



(1.5) 
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However, we get many other (/-analogues of r F s by adding upper or lower parameters to 
the left hand side of (1.5) which are equal to or which depend on q in such a way that 
they tend to a limit ^ 1 as q f 1. Note that the notation (1.3) has the drawback that some 
rescaling is necessary before we can take limits for q — > 1. On the other hand, parameters 
can be put equal to zero without problem in (1.3), which would not be the case if we 
worked with a±, . . . , a r , bi, . . . , b s as in (1.5). In general, (/-hypergeometric series can be 
studied in their own right, without much consideration for the q = 1 limit case. This 
philosophy is often (but not always) reflected in the notation generally in use. 

It is well known that the confluent \F\ hypergeometric function can be obtained from 
the Gaussian 2 F\ hypergeometric function by a limit process called confluence. A similar 
phenomenon occurs for (/-series. Formally, by taking termwise limits we have 



lim 



a-i, 



' °, r ; q, — 

, S CL r 



ai, . . . , a r - 

&1,... ,bg 



(1.6) 



In particular, this explains the particular choice of the factor ((— l) k (/ fc ( fc_1 )/ 2 ) 1+s r in 
(1.3). If r = s + 1 this factor is lacking and for r < s + 1 it is naturally obtained by the 
confluence process. The proof of (1.6) is by the following lemma. 

Lemma 1.1 Let the complex numbers a^, k G Z + , satisfy the estimate | ctfe j < R~ k for 
some R > 0. Let 



F(b;q,z) := Va fc (b;q) k z k , 



fc=0 



\z\ < R, 6gC, 0<g<l. 



Then 



lim F(b; q, z/b) = V a k (-l) k q k ^-^/ 2 z k , 

h — vnn ' * 



(1.7) 



k=0 

uniformly for z in compact subsets of C. 

Proof For the kth term of the series on the left hand side of (1.7) we have 

c fc :=a k (6; q)k {z/b) k 

=a k (b- 1 - 1) (b- 1 -q)... (b- 1 - q k ~ l ) z k . 
Now let \b\ > 1 and let E N be such that q 1 *' 1 > > q N . If k < N then 

\c k \ < \a k \q k ^ k - 1 ^ 2 {2\z\) k . 



If k > N then 



Kb- 1 - 1X&" 1 -q)... (b~ l - q k ~ l )\ < 2 k q N ( N ~V/ 2 \b\ N ~ k 

< 2 fc| 6 |(l-iV)/2| 6 |iV-fc 

< 2 k \b\~ k/2 , 



so 



|c fc |<|a fc |(2|^||6|- 1 / 2 ) fc . 
Now fix M > 0. Then, for each e > we can find K > and B > 1 such that 



| c fe | < e if \z\ < M and |6| > B. 



k=K 



Combination with the termwise limit result completes the proof. 



□ 
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It can be observed quite often in literature on (/-special functions that no rigorous 
limit proofs are given, but only formal limit statements like (1.6) and (1.5). Sometimes, 
in heuristic reasonings, this is acceptable and productive for quickly finding new results. 
But in general, I would say that rigorous limit proofs have to be added. 

Any terminating power series 



E 

k=0 



c k z 



can also be written as 



^c n _ fc (1/z)' 



k=0 



When we want to do this for a terminating (/-hypergeometric series, we have to use that 



(a;q) n -k = 



(a; q) r 



(q^a-q-^k 



= {a;q)n 



(_!)fcgfc(fc-i)/2 ( a -y-ny 
(q 1 - n a- 1 ;q) k 



and 



(g n ;q)n-k = (q n ;q) n (q n ;q h _ ( _-,\ n - n (n+i)/2 (g k (n+l)k 



(q;q) n -k (?;?)n (q 1 ;q x )k 

Thus we obtain, for n G Z+, 



S+10S 



a±, . . . , a s 



6i,...,6 a 



;q,z 



_ (_l) n q- n ( n + 1 )/ 2 ( fl l7 • • • ; a s'i q)n 



q , \q- 



-n+i„-i 



(b 1 ,...,b s ;q) n 

1 b^\...,q- n+1 b- 1 _ ? n+1 6i...6 s 
. . , q-'^^a b 



n > 



-n+l r ~l ' ,q > 



ai . . . a s z 



(1.8) 



Similar identities can be derived for other r cj) s and for cases where some of the parameters 
are 0. 

Thus, any explicit evaluation of a terminating (/-hypergeometric series immediately 
implies a second one by inversion of the direction of summation in the series, while any 
identity between two terminating (/-hypergeometric series implies three other identities. 



1.2. The (/-integral. Standard operations of classical analysis like differentiation and 
integration do not fit very well with (/-hypergeometric series and can be better replaced by 
(/-derivative and (/-integral. 

The q-derivative D q f of a function / on an open real interval is given by 



(1.9) 



and (D q f)(0) := f'(0) by continuity, provided /'(0) exists. Note that lim q ^i(D q f)(x) 
= f'(x) if / is differentiable. Note also that, analogous to d/dx (1 — x) n = —n (1 — x) n_1 , 
we have 

1 - Q n 

f(x) = Or; q) n (£>,/)(*) = "y^y (^; q)n-i- (1-10) 
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Now recall that < q < 1. If D q F = f and / is continuous then, for real a, 

oo 

F(a) - F(0) =a(l-q)J2 /(«?*) 9*- (1-H) 



A;=0 

This suggests the definition of the q-integral 

pa 00 

/ /(*) d,* := a (1 - g) ^T/(ag fc )g fc . (1.12) 

Note that it can be viewed as an infinite Riemann sum with nonequidistant mesh widths. 
In the limit, as q f 1, the right hand side of (1.12) will tend to the classical integral 
Jo f(x)dx. 

From (1.11) we can also obtain F(a) — F(b) expressed in terms of /. This suggests 
the definition 

pb pa pb 

\ f(x)d q x:= / f(x)d q x- / f(x)d q x. (1.13) 

Jo JO JO 

Note that (1.12) and (1.13) remain valid if a or b is negative. 

There is no unique canonical choice for the g-integral from to oo. We will put 



poo 00 

/ f{x) d q x := (1 - q) f^)^ 

Jo u - 



k= — oo 

(provided the sum converges absolutely). The other natural choices are then expressed by 

poo 00 

a / f(ax)d q x = a(l-q) V f(aq k )q k , a > 0. 

J ° k=-oo 

Note that the above expression remains invariant when we replace a by aq n (n 6 Z). 
As an example consider 

«1 oo . 

/ x a d q x=(l-q)J2<l Ha+1) = i Ik- Re«>-1, (1.14) 

which tends, for q j 1, to 

1 r 1 

x Q &. (1.15) 



CK+ 1 







From the point of view of explicitly computing definite integrals by Riemann sum ap- 
proximation this is much more efficient than the approximation with equidistant mesh 
widths 



-±{- 



n * — ' \ n 

k=l 



The ^-approximation (1.14) of the definite integral (1.15) (viewed as an area) goes essen- 
tially back to Fermat. 
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1.3. Elementary examples. The q-binomial series is defined by 

^o(a;-;^):=f;^A \z\ < 1. (1.16) 

Here the '— ' in a r s expression denotes an empty parameter list. The name "(/-binomial" 
is justified since (1.16), with a replaced by q a , formally tends, as q | 1, to the binomial 
series 

1 Fo(a;^):=f;^^ = (l-^)- a , |z| < 1. (1.17) 



The limit transition is made rigorous in [21, Appendix A]. It becomes elementary in the 
case of a terminating series (a := q~ n in (1.16) and a := — n in (1.17), where n G Z+). 
The (/-analogue of the series evaluation in (1.17) is as follows. 



Proposition 1.2 We have 

^(n- -■ n -y\ — — _ . 



1 Ma;-,q,z)= i -^^, \z\<l. (1.18) 



In particular, 

1 Mq- n ---q,z) = {q- n z;q) n . (1.19) 
Proof Put /i a (z) := i0o(a; — ; z). Then 

(1 -z)/i (z) = (1 -az)h a (qz), hence ft a (z) = -y- — h a (qz). 

(1 - ^) 



Iteration gives 



(*;?)n 



Now use that h a is analytic and therefore continuous at and use that h a (0) = 1. Thus, 
for n — > oo we obtain (1.18). □ 

The two most elementary (/-hypergeometric series are the two q-exponential series 

°° z k 1 

e q (z) := !0 O (O; -; «,*) = £ 7^T = JZT^i \ z \ < X > (1.20) 

^ {q,q)k {z,q)oo 

and 

°^ q k(k-i)/2 z k 

E q (z) := (f) (-;-;q,-z) = 2_^ — r = (- Z ;q) 00 , z G C. (1.21) 

fc=o [q,q)k 

The evaluation in (1.20) is a specialization of (1.18). On the other hand, (1.21) is a 
confluent limit of (1.18): 

lim i0 o (a; -; 9, -z/a) = o 0o(-; -; 9, -^) 
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(cf. (1.5)). The limit is uniform for z in compact subsets of C, see Lemma 1.1. Thus the 
evaluation in (1.21) follows also from (1.18). 
It follows from (1.20) and (1.21) that 

e q (z)E q (-z) = 1. (1.22) 

This identity is a (/-analogue of e z e~ z = 1. Indeed, the two (/-exponential series are 
(/-analogues of the exponential series by the limit formulas 

\imE q ((l — q)z) = e z = lime g ((l — q)z). 

<7tl <7T1 

The first limit is uniform on compacta of C by the majorization 

qk{k-i)/2 _ ~ k 



< 



k\ 



The second limit then follows by use of (1.22). 

Although we assumed the convention < g < 1, it is often useful to find out for 
a given (/-hypergeometric series what will be obtained by changing q into q~ x and then 
rewriting things again in base q. This will establish a kind of duality for (/-hypergeometric 
series. For instance, we have 

e q -i{z) = E q (-qz), 

which can be seen from the power series definitions. 

The following four identities, including (1.20), are obtained from each other by trivial 
rewriting. 



E 



(1 - q) k z k 1 



k _ A] (q;q)k {{i- q)z- 1 q) 00 " 



j, 



x ^ z K 1 

^ Q (q;q)k (z;q)c 



oo 



fl -n) 1 - b a kb (a k+1 -a) = 

( q) h {q ,q)o ° (^)oo(i-^- 1 ' 

^((1-^;^^=^-^^, Re6>0. (1.23) 



As q | 1, the first identity tends to 

OO h. 

~ rv 

k=0 

while the left hand side of (1.23) tends formally to 
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Since this last integral can be evaluated as T(b), it is tempting to consider the right hand 
side of (1.23) as a the q-gamma function. Thus we put 



(q; V)c 



r ' (z):= (^r^ = X ^(-d-^XM. «»><>■ 

where the last identity follows from (1.23). It was proved in [21, Appendix B] that 

\imT q (z) =T(z), z ^ 0,-1,-2,... . (1.24) 

We have just seen an example how an identity for (/-hypergeometric series can have 
two completely different limit cases as q | 1. Of course, this is achieved by different 
rescaling. In particular, reconsideration of a power series as a (/-integral is often helpful 
and suggestive for obtaining distinct limits. 

Regarding T q it can yet be remarked that it satisfies the functional equation 

r,(* + i) = ^r,(z) 

and that 

r> + i) = i^, nez + . 

Similarly to the chain of equivalent identities including (1.20) we have a chain including 
(1.18): 

g (q a ;q) k z k _ (q a z;q) 00 



(q;q) k (::(/) x 

(1 - a) V o kb = (1-9) ^t + \q) 

^ q ((/ fc+a ;(/)oc (q a ,q b ;q)c 



k=0 



7 (g a t; ? )oo q T q (a + b) K J 



In the limit, as q f 1, the first identity tends to 



fc=0 



while (1.25) formally tends to 



1 t b - 1 (l-t) a - 1 dt=^^ • (1-26) 
Thus (1.25) can be considered as a (/-beta integral and we define the q-beta function by 



„,„ M . r .(") r »C') ('-«) to.g° +t ;9)oo /Vi fat; 9)00 , , 97 > 
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1.4. Heine's 2 <fri series. Euler's integral representation for the 2 Fi hypergeometric 
function 

= m$ZT) J! ^ (1 ~ ^ P " ^ *■ (1.28) 
Re c> Re 6 > 0, | arg(l - z) | < 7r, 

(cf. [13, 2.1(10)]) has the following (/-analogue due to Heine. 

, (q b c x r <?( c ) A 1 +6-1 (*?;9)oo (tzq a ;q) QO 

'* (, '' ;,;l2)= r,( W -l)i, ' (*«-»;«)„ ~5^r **•*" (1-29) 

Re6 > 0, |z| < 1. 

Note that the left hand side and right hand side of (1.29) tend formally to the corre- 
sponding sides of (1.28). The proof of (1.29) is also analogous to the proof of (1.28). 
Expand (tzq a ;q) 00 /(tz;q) 00 as a power series in tz by (1.16), interchange summation and 
(/-integration, and evaluate the resulting (/-integrals by (1.25). 

If we rewrite the (/-integral in (1.29) as a series according to the definition (1.12), and 
if we replace (/ a , q b , q c by a, 6, c then we obtain the following transformation formula: 

2^1 (a, b; c; q, z) = ^ Z,q ^°° h 2 (f>i(c/b, z; az; q, b). (1.30) 

Although (1.29) and (1.30) are equivalent, they look quite different. In fact, in its form 
(1.30) the identity has no classical analogue. We see a new phenomenon, not occurring for 
2-F1, namely that the independent variable z of the left hand side mixes on the right hand 
side with the parameters. So, rather than having a function of z with parameters a, 6, c, 
we deal with a function of a, 6, c, z satisfying certain symmetries. 

Just as in the classical case (cf. [13, 2.1(14)]), substitution of some special value of z in 
the (/-integral representation (1.29) reduces it to a (/-beta integral which can be explicitly 
evaluated. We obtain 

20i (a, ft; c; g ,c/(a&))= , |c/(a6)| < 1, (1.31) 

(c, c/(ab); q)^ 

where the more relaxed bounds on the parameters are obtained by analytic continuation. 
The terminating case of (1.31) is 

2 MQ-\b-c-q,cq n /b)= { -^^, n G Z+. (1.32) 

l C 5 q)n 

The two fundamental transformation formulas 



2 Fi(a, b- c; z) = (l- z)~ a 2 F 1 (a, c-b;c; z/(z - 1)) 
= (l-z) c - a - b 2 F 1 (c-a,c-b; c; z) 



(1.33) 
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(cf. [13, 2.1(22) and (23)]) have the following (/-analogues. 



2 (f)i(a,b;c;q,z) 



(az;q)c 

(abz/c; q) 



2^2(0, c/b; c, az; q, bz) 

20i (c/a, c/b; c; q, abz/c). 



(1.34) 
(1.35) 



Formula (1.35) can be proved either by threefold iteration of (1.30) or by twofold iteration 
of (1.34). The proof of (1.34) is more involved. Write both sides as power series in z. Then 
make both sides into double series by substituting for (b; q)^/ (c; q)k on the left hand side a 
terminating 20i (cf. (1-32)) and by substituting for (aq h z; q) OQ /(z; q)oo on the right hand 
side a (/-binomial series (cf. (1.16)). The result follows by some rearrangement of series. 
See [15, §1.5] for the details. 

Observe the difference between (1.33) and its (/-analogue (1.34). The argument z/(z — 
1) in (1.33) no longer occurs in (1.34) as a rational function of z, but the ^-variable is 
distributed over the argument and one of the lower parameters of the 202 • Also we do not 
stay within the realm of 20i functions. 

Equation (1.8) for s : = 1 becomes 



(b;q) n . ,„ ± (__ n q -n+l_ q -n+l^ q n+l c 



2 0! (q~\ 6; c; q, z) = q -<*«+W ^ (-*)» 2 fc ( q 



c b 



bz 



n G Z+. 

(1.36) 

We may apply (1.36) to the preceding evaluation and transformation formulas for 20i in 
order to obtain new ones in the terminating case. From (1.32) we obtain 



, . _ n , (c/b;q) n b n 

20i [q ,6;c;g,g) = — r , n E Z+. 

(c; q)n 



(1.37) 



By inversion of direction of summation on both sides of (1.34) (with a := q n ) we obtain 



/ / —n 1 \ (c b;q) n 

201 [q , b; c; q, z) = — r — 3 02 

c; q)n 



q n ,b,bzq n /c 
bq 1 ~ n /c,0 



n e Z+. 



(1.38) 



A terminating 2 0i can also be transformed into a terminating 3 2 with one of the 
upper parameters zero (result of Jackson): 



20i (q n , b; c; <?, z) = (q n bz/c; q) n 3 2 



Q" n ,c/6,0 
c, cqb~ 1 z~ 1 



(1.39) 



This formula can be proved by applying (1.19) to the factor (cq k+1 b~ 1 z~ 1 ; q) n -k occuring 
in the kth. term of the right hand side. Then interchange summation in the resulting double 
sum and substitute (1.37) for the inner sum. See [20, p. 101] for the details. 
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1.5. A three-term transformation formula. Formula (1.36) has the following gen- 
eralization for non-terminating 2<Pi'- 

(a,q/c,c/b,bz/q,q 2 /bz;q) oc ( aq bq q 2 

2<Pi{a,b;c;q,z) + y—. JTl—/ 771 S \ 201 — , — ; — ;q,z 

_ (abz/^q/^aq/^cq/jabz)^)^ ^ ( aq aq ^ eg \ cq < < l 
{bz/c,q/b,aq/c,cq/{bz)]q)oo \ ' c' 6' 7 abzJ ' ab 

This identity is a (/-analogue of [13, 2.1(17)]. In the following proposition we rewrite (1.40) 
in an equivalent form and next sketch the elegant proof due to [27]. 

Proposition 1.3 Suppose {q n | n G Z + } is disjoint from {a~ 1 q~ n , b~ 1 q~ n | n G Z + }. 
Then 

(z, gz" 1 ; 2 fc(a, 6; c; 9 , z) = (az, (/ert" 1 ; g) (c/a ' b] qU 



(c,6/a; g)oo 

X20i(a, (/a/c; (/a/6; (/, (/c/ (a&z)) + (a < — ► b). (1.41) 



Proof Consider the function 

(a, 6, cw,q, qwz~ x ^w^z^q)^ 1 



F(w) 



(aw,bw,c,w 1 ;q)oo w 



Its residue at q n (n G Z + ) equals the n th term of the series on the left hand side of (1.41). 
The negative of its residue at a~ 1 q~ n (n G Z + ) equals the n th term of the first series on 
the right hand side of (1.41), and the residue at b~ 1 q~ n is similarly related to the second 
series on the right hand side. Let C be a positively oriented closed curve around in C 
which separates the two sets mentioned in the Proposition. Then (27ri) _1 f c F(w) dw can 
be expressed in two ways as an infinite sum of residues: either by letting the contour shrink 
to {0} or by blowing it up to {oo}. EH 

If we put z := q in (1.40) and substitute (1.31) then we obtain a generalization of 
(1.37) for non-terminating 2</>i: 

, fh \, (a,M/c;g)oo f aq bq q 2 \ (abq/^q/c^q)^ 
201 (o, b; c; q, q) + — — — 2 0i — , — ; — 5 9, 9 = 7 / Z / r - • ( L42 ) 



(aq/c,bq/c,c/q;q)oo \c c c J (aq/c,bq/c;q) 
By (1.13) this identity (1.42) can be equivalently written in (/-integral form: 

(ct, qt; q)oq d t _ (1 - q) (abq/c, q/c, c, q; q)^ ^ ^ 



qc 



(at,bt',Q)oo (aq/c 7 bq/c,a 7 b;q) c 



Replace in (1.43) c, a, b by — (/ c , —9", (/ b , respectively, and let q | 1. Then we obtain 
formally 



+ «)•-■ (i ^7 it 1 ^ - (L44) 

.1 1 (a + — c + 1) 

Note that, although it is trivial to obtain (1.44) from (1.26) by an affine transformation 
of the integration variable, their (/-analogues (1.25) and (1.43) are by no means trivially 
equivalent. 
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1.6. Bilateral series. Put a := q in (1.40) and substitute (1.16). Then we can combine 
the two series from to oo into a series from — oo to oo without "discontinuity" in the 
summand: 



(b;q)k k (q,c/b,bz,q/(bz);q) c 

2^ - — ~ z - 



^(c;q)k (c,q/b,z,c/(bz);q) 
Here we have extended the definition of (6; q)k by 

(M)oo 



\c/b\ < \z\ < 1. 



(1.45) 



(6; o)k 



(bq k ;q) c 



Formula (1.45) was first obtained by Ramanujan (Ramanujan's iV'i -summation formula). 
It reduces to the (/-binomial formula (1.18) for c := q n (n G Z + ). In fact, this observation 
can be used for a proof of (1.45), cf. [17], [2, Appendix C]. Formula (1.45) is a (/-analogue 
of the explicit Fourier series evaluation 

E (ft)fc jke _ r(c) r(i - 6) i(1 _ c)(:e ,_ 7r ) _ i0s C _ b _i 



k= — oo 



r( c - 6) 



< 9 < 2tt, Re(c-6- 1) > 0, 



where (b) k := T(b + k)/T(b), also for k = -1, -2, 
Define bilateral q-hypergeometric series by 



ai, . . . , a r 



b 1 ,...,b l 



E (ai, . . . , a r ; g)fc , fc k(k-i)/2\ s - r k 
(h,...,b s ;q) k (i 1} q 



k= — oo 



~-r+l(Ps 



(1.46) 



ai, . 



+ 



(&i - g) ■ ■ ■ (b s - q) 
(ai - 5) . . . (a r - g) z 



q 2 /b u ...,q 2 /b s ,q 



h...b s 



_q 2 /ai, . . . , q 2 /a r , 0, . . . , 0' ai . . . a r z 
where ai, . . . a r , b\, . . . , b s ^ and s > r. The Laurent series in (1.46) is convergent for 

h...b 



(1.47) 



and for 



<i\ . . . a r 
b x ...b 



<\z\ if s > r 



a\ . . .a r 



< \z\ < 1 if s 



If some of the lower parameters in the r ip s are then a suitable confluent limit has to be 
taken in the s +i0 s of (1.47). 

Thus we can write (1.45) as 



iV>i(&;c;g,z) = 



(g,c/6, bz,q/(bz);q) c 
{c,q/b,z,c/(bz)\q) (> 



\c/b\ < \z\ < 1. 



(1.48) 
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Replace in (1.48) z by z/b, substitute (1.47), let b — > oo, and apply (1.6) Then we obtain 

qk(k-i)/2 z k (q jZ ,q/z;q) c 

k= — oo 



o^i(-; c; q, z) := ^ 



\z\ > c 



(1.49) 



(c;g)fc (c,c/^; q)oo 

In particular, for c = 0, we get the Jacobi triple product identity 

oo 

oVi(-;0; E (-l) fc ^ (fc - 1)/2 ^ = (9,^,^;9)oc, z^O. (1.50) 

fc= — oo 

The series in (1.49) is essentially a theta function. With the notation [14, 13.19(9)] 
we get 

oo 

6 4 (x;q) := ^ (-l) fe g fc2 e 2 ™ fc:B 

k= — oo 

=oVi(-;0;(/ 2 ,(/e 2 ™) 



= J] (! - 4 2fc ) (1 - V"" 1 cos(2ttx) + </ 4fc - 2 ), 
fe=i 

and similarly for the other theta functions 6i(x;q) (i = 1, 2, 3). 

1.7. The (/-hypergeometric (/-difference equation. Just as the hypergeometric dif- 
ferential equation 

z(l-z) u"{z) + (c-(a + b + l)z) u'{z) - abu(z) = 
(cf. [13, 2.1(1)]) has particular solutions 

ui(z) := 2-Pi(a, 6; c; z), u 2 (z) := z 1 ~ c 2-^1(0 — c + 1, 6 — c + 1; 2 — c; z), 
the q-hyper geometric q-difference equation 



z(q c 



„a+b+l. 



1 - q c 
1-9 



6+1 



l-q 



1-q 

l-q a l-q b 



(D q u)(z) 
u(z) = (1.51) 



has particular solutions 



iti(z) := 2 (f>i(q a ,q b ;q c ;q,z), 

-u 2 (^j:=2 20i {q ,q ,q ,q,z). 



(1.52) 
(1.53) 



There is an underlying theory of (/-difference equations with regular singularities, similarly 
to the theory of differential equations with regular singularities discussed for instance in 
Olver [29, Ch. 5]. It is not difficult to prove the following proposition. 
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Proposition 1.4 Let A(z) := J2T=o ak zk an< ^ B(z) '■— Y^k=o^kZ k be convergent power 
series. Let A G C be such that 

(1 - q X+k ) (1 - q^- 1 ) l~q X+k f = 0, k = 0, , . 

+ a o — i ; — ^ b o s , n l _ i o U- b4 J 



(1-g) 2 U 1-9 U 1^0, fc = l,2,.... 

Then the (/-difference equation 

{D 2 q u){z) +zA(z) (D q u)(z) + B(z)u(z) = (1.55) 
has an (up to a constant factor) unique solution of the form 



(z) = £>z A+ *. (1.56) 



fc=0 



Note that (1.51) can be rewritten in the form (1.55) with a = (q~ c — 1)/(1 — g), 
bo = 0, so (1.54) has solutions A = and — c + 1 (mod (2ni logg _1 )Z) provided c ^ Z 
(mod (2/u logg _1 )Z). For the coefficients c k in (1.56) we find the recursion 



c k+1 (1 - q a + x + k ) (1 - q b+W 



C k (l-q c + X + k ) (1-gA+fc+l) 



Thus we obtain solutions iti, -U2 as given in (1.52), (1.53). 

Proposition 1.4 can also be applied to the case z = oo of (1.51). Just make the 
transformation z h- > z _1 in (1.51). One solution then obtained is 



w 3 (z) := ^ a 20i 
Now (1.40) can be rewritten as 



rl a r ,a— c+1 

y ' y . „ „-a-b+ c +i -i 
-&+i >g>g 



(g^V c ,g c ^g)^ (g b ^,g 2 b z ^9)00^ 1 , . 

Ul[Z) (y c - 1 ,y a - c+1 ,y 1 - b ;y)oo (g^g-^" 1 ;^ WJ 

= (g^g^+^gjoo (g 0+b - c z,g c - a - b+1 z- 1 ;g)ocz° , , , , 

(? l-6 )9 a-c + l ;?)oo ( g6-c Zj gC-6+l z -l ;?)oo ^j- U^j 

Note that 113 is not a linear combination of iti and 112, as the coefficients of and W3 in 
(1.57) depend on z. However, since an expression of the form 

^ (g a z,g 1 - a z- 1 ;g) 00 z a - /3 
(qPz, g 1_/3 2: _1 ;g)oo 

is invariant under transformations 2; 1— > qz, each term in (1.57) is a solution of (1.51). Thus 
everything works fine when we restrict ourselves to a subset of the form {zq q k \ k G Z}. 

The g-analogue of the regular singularity at 2 = 1 for the ordinary hypergeometric 
differential equation has to be treated in a different way. We can rewrite (1.51) as 

(q c - q a+b z) u(qz) + (-g c - q + (q a + q h )z) u(z) + (q - z) u{q- x z) = 0. 
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It can be expected that the points z = q c ~ a ~ b and z = q, where the coefficient of u(qz) 
respectively u(q~ 1 z) vanishes, will replace the classical regular singularity z = 1, see also 
(1.31), (1.37) and (1.42). One can find solutions related to these singularities, for instance 



U 4 (z) := 302 
U 5 (z) := Z~ b 302 



n a n b n a+b—c~ 

q i q i q 



ga+b-c+l^Q 

q b , q b ~ c+ \ qz 



q a+b - c+1 ,0 

and many others. A systematic theory for such singularities has not yet been developed. 
1.8. (/-Bessel functions. The classical Bessel function is defined by 



00 (-1) k 

k=o v ; 



cf. [14, Ch. 7]). Jackson (1905) introduced two g-analogues of the Bessel function: 

q) := ( ^ 1 \ 9) °° (z/2)" 2 0i(O, 0; <f + 1 ; q, -z 2 /4) 
\Qi Q)oo 

=(z/2) 



(1.58) 



4 2 \z; q) := ( ^ 1 \ 9) °° (z/2)" *i(-; «, V/4) 



loo 
oo 



/ /2) „ (^ +fc+1 ;9)oo 9 fc(fc+I/) (~i) fc (^/2) 2fc 
fr («;«)«, 



(1.59) 



Formally we have 



KmJ<?\(l-q)z;q) = J v (z), i = 1, 2 



(cf. (1.24)). The two (/-Bessel functions jj> (z; q) can be simply expressed in terms of each 
other. From (1.35) and (1.6) we obtain 

2 0i(O,6;c;g,z) 



A further confluence with 6^0 yields 



(z;q)c 



1 



^(c/b;c;q,bz). 



20i (0, 0; c; (/, z) = — o 0i(-; c; g, cz). 
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Formally we have 



This can be rewritten as 

4 2 Hz; q ) = (-z 2 /4; q U4 1 \z; q ). 
Yet another (/-analogue of the Bessel function is as follows. 

Mz; q) := (g ^ 1 \ g) °° z v ^(0; ? , ^ 2 ) 

~ C^+Hl. „\ (_-|\fc „fc(fc+l)/2 „2fc 

= {Q , ( lj q r — . (1.60) 

f-^ (g; 9)00 (g;g)fc 

limJ I/ ((l-g 1 / 2 )z;g) = J„(z). 

This (/-Bessel function is not simply related to the (/-Bessel functions (1.58), (1.59). As far 
as we know, it was first introduced by Hahn (in a special case) and by Exton (in general). 
In recent work by Koornwinder & Swarttouw [23] a satisfactory (/-analogue of the Hankel 
transform could be given in terms of the (/-Bessel function (1.60). 

1.9. Various results. Goursat's list of quadratic transformations for Gaussian hyper- 
geometric functions can be found in [13, §2.11]. In a recent paper Rahman & Verma [30] 
have given a full list of (/-analogues of Goursat's table. However, all their formulas involve 
on at least one of both sides an 8</>7 series. Moreover, a good foundation from the theory of 
(/-difference equations is not yet available. For terminating series many of their §^7 's will 
simplify to 4^3 's. In section 2 we will meet some natural examples of these transformations 
coming from orthogonal polynomials. 

An important part of the book by Gasper & Rahman [15] deals with the derivation of 
summation and transformation formulas of s +i<^ s functions with s > 1. A simple example 
is the q-Saalschiitz formula 

3 Ma,b,q- n -,c,abc-V- n -, q , q )= ^ c/a ;^ q)n , neZ + , (1.61) 

(c, c/{ab);q) n 

which follows easily from (1.35) by expanding the quotient on the right hand side of (1.35) 
with the aid of (1.18), and next comparing equal powers of z at both sides of (1.35). 
Formula (1.61) is the (/-analogue of the Pfaff-Saalschiitz formula 

J-! / 7 -1 1 . , (c- a) n (c-b) n 
3 F 2 (a, b, -n; c, 1 + a + b - c - n;l) = — — ; — , n G Z + , 

(c)„ (C - a - b) n 

which can be proved in an analogous way as (1.61). 

An example of a much more involved transformation formula, which has many impor- 
tant special cases, is Watson's transformation formula 



a, qa 



1 / 2 ,- ? a 1 / 2 ,6,c,d,e, 9 - n aV +n 



; q 

a l / 2 ,—a l / 2 ,aq/b,aq/c,aq/d,aq/e,aq n+1 ' ' bcde 



(ag, aq/(de);q) n 
(aq/d, aq/e;q) n 



q n ,d,e,a<?/(6c) 
aqjb.aqj c,deq~ n I a ' 



neZ + , (1.62) 
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cf. Gasper & Rahman [15, §2.5]. Then, for a 2 q n+1 = bcde, the right hand side can be 
evaluated by use of (1.61). The resulting evaluation 

(ag, ag/(6c), aq/(bd), aq/(cd); q) n 



(aq/b, aqj c, aqj d, aqj (bed) ; q) n 



of the left hand side of (1.62) subject to the given relation between a, 6, c, <i, e, q is called 
Jackson's summation formula, cf. [15, §2.6]. 
The famous Roger s-Ramanuj an identities 

, , ft x f q e (g 2 ,g 3 ,g 5 ;g 5 )oo 

o 0i(-;O;g,g := > — = r , 

g;g -t g;g u- 



fc=0 



, , n 2 , ^ g fc(fc+1) (g,g 4 ,g 5 ;g 5 )c 

o0i(-;O;g,g ) := ^ -— - = — — 



,_ (g;g)fc (g;g) . 

have been proved in many different ways (cf. Andrews [2]), not only analytically but also 
by an interpretation in combinatorics or in the framework of Kac- Moody algebras. A quick 
analytic proof starting from (1.62) is described in [15, §2.7]. 

Exercises to §1 

1.1 Prove that 

(a;g) n = (-arg^" 1 )/ 2 (a^g 1 -; q) n . 

1.2 Prove that 

(a;g) 2 n = (a; g 2 ) n (ag;g 2 ) n , 
(a 2 ;g 2 ) n = (a; g) n (-a;g) n , 
(a; q) OQ = (a 1 ' 2 , -a 1 / 2 , (ag) 1 ^, -(ag) 1 ^ ?)<x) . 

1.3 Prove the following identity of Euler: 

(-g;g)oo (g;g 2 )oo = 1. 

1.4 Prove that 

~ /_ 1 y-m ? (Z-m)(Z-m-l)/2 

, ^ r„(n - < + 1) r„(< - m + 1) 

Do it first for g = 1. Start for instance with e z e~ z = 1 and, in the g-case, with 
e q (z)E q (-z) = l. 
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1.5 Let 

^ oo 

7 — r~ = y2 a nQ n 

fro 

be the power series expansion of the left hand side in terms of q. Show that a n equals the 

number of partitions of n. 

Show also that the coefficient b k , n in 



k 00 



is the number of partitions of n with highest part k. Give now a partition theoretic proof 
of the identity 



(9; 9)oo t^iVirik 

1.6 In the same way, give a partition theoretic proof of the identity 

1 00 „mfc 

7 v— = / \ , men. 

(? m ;?)oo £~' (vq)k 

1.7 Give also a partition theoretic proof of 

00 fc(fc-l)/2 mfc 

fro 

(Consider the problem first for m = 1.) 

1.8 Let GF(p) be the finite field with p elements. Let A be a n x n matrix with entries 
chosen independently and at random from GF(p), with equal probability for the field 
elements to be chosen. Let q := p~ x . Prove that the probability that A is invertible is 
(q;q) n - (See SIAM News 23 (1990) no.6, p.8.) 

1.9 Let GF(p) be as in the previous exercise. Let V be an n-dimensional vector space 
over GF(p). Prove that the number of /c-dimensional linear subspaces of V equals the 
(/-binomial coefficient 



n 



(P',P)n 



kip' (p;p) k (p;p)n-k' 
1.10 Show that 

n 

(ab;q) n = ^2^] b k (a;q) k (M)n-fc- 
fc=o q 

Show that both Newton's binomial formula and the formula 



(a + b) n = ^2 ( \ ) ( a )k K-t 



fe=0 

are limit cases of the above formula. 
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2. (/-Analogues of the classical orthogonal polynomials 

Originally by classical orthogonal polynomials were only meant the three families of Ja- 
cobi, Laguerre and Her mite polynomials, but recent insights consider a much bigger class 
of polynomials as "classical". On the one hand, there is an extension to hypergeometric 
orthogonal polynomials up to the 4F3 level and including certain discrete orthogonal poly- 
nomials. These are brought together in the Askey tableau, cf. Table 1. On the other hand 
there are (/-analogues of all the families in the Askey tableau, often several g-analogues 
for one classical family (cf. Table 2 for some of them). The master class of all these q- 
analogues is formed by the celebrated Askey- Wilson polynomials. They contain all other 
families described in this chapter as special cases or limit cases. Good references for this 
chapter are Andrews & Askey [4] and Askey & Wilson [8]. See also Atakishiyev, Rahman 
& Suslov [10] for a somewhat different approach. 

Some parts of this section contain surveys without many proofs. However, the subsec- 
tions 2.3 and 2.4 on big and little (/-Jacobi polynomials and 2.5 and 2.6 on the Askey- Wilson 
integral and related polynomials are rather self-contained. 

2.1. Very classical orthogonal polynomials. An introduction to the traditional 
classical orthogonal polynomials can be found, for instance, in [14, Ch. 10]. One possible 
characterization is as systems of orthogonal polynomials {p n }n=o,i,2,... which are eigenfunc- 
tions of a second order differential operator not involving n with eigenvalues A n depending 
on n: 

a(x)p'^(x) + b(x)p' n (x) + c(x)p n (x) = X n p n (x). (2.1) 

Because we will extend the definition of classical orthogonal polynomials in this section, we 
will call the orthogonal polynomials satisfying (2.1) very classical orthogonal polynomials. 
The classification shows that, up to an affine transformation of the independent variable, 
the only cases are as follows. 

Jacobi polynomials 

P^\x) := ^ a+ ^ n 2 F 1 (-n,n + a + p+ 1; a +l;(l-x) /2), a,p>-l, (2.2) 
n! 

orthogonal on [—1, 1] with respect to the measure (1 — x) a (1 + x)@ dx; 
Laguerre polynomials 

L%(x) := ^±i^ 1 F 1 (-ri;a + l;a;), a > -1, (2.3) 

orthogonal on [0, 00) with respect to the measure x a e~ x dx; 
Hermite polynomials 

H n (x) := (2x) n 2 F (-n/2, (1 - n)/2; -; -x" 2 ), (2.4) 

2 

orthogonal on (—00, 00) with respect to the measure e~ x dx. 
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In fact, Jacobi polynomials are the generic case here, while the other two classes are 
limit cases of Jacobi polynomials: 

L$(x) = lim P^\l-2x/{3), (2.5) 
H n (x) = 2 n n\ lim a~ n / 2 P^ia^^x). (2.6) 

Hermite polynomials are also limit cases of Laguerre polynomials: 

H n {x) = (-l) n 2 n l 2 nl lim a~ n/2 L«((2a) 1/2 x + a). (2.7) 

a^oo 

The limit (2.5) is immediate from (2.2) and (2.3). The limit (2.6) follows from (2.4) and a 
similar series representation for Gegenbauer or ultraspherical polynomials (special Jacobi 
polynomials with a = (3): 

P^\x) = {a + ( l )n ^ l [ 2)n (2x) n 2 Fi(-n/2, (1 - n)/2; -a - n + 1/2; x~% 
{la + l) n n\ 

cf. [14, 10.9(4) and (18)]. The limit (2.7) cannot be easily derived by comparison of series 
representations. One method of proof is to rewrite the three term recurrence relation for 
Laguerre polynomials 

(n + 1) L% +1 (x) - (2n + a + 1 - x) L™{x) + (n + a) L^x) = 

(cf. [14, 10.12(8)]) in terms of the polynomials in x given by the right hand side of (2.7) 
and, next, to compare it with the three term recurrence relation for Hermite polynomials 
(cf. [14, 10.13(10)]) 

H n+1 (x) -2xH n {x) + 2nH n _ 1 (x) =0. 

Very classical orthogonal polynomials have other characterizations, for instance by 
the existence of a Rodrigues type formula or by the fact that the first derivatives again 
form a system of orthogonal polynomials (cf. [14, §10.6]). Here we want to point out that, 
associated with the last two characterizations, there is a pair of differential recurrence 
relations from which many of the basic properties of the polynomials can be easily derived. 
For instance, for Jacobi polynomials we have the pair 

i = n+a ^ +1 et wi) w, (2.8) 

(1 - x)-" (1 + x)-' 3 ^- ((1 - x) a+1 (1 + x f +1 P^°t 1,,3+I, (a-)) = -2nPi"' 0> (x). (2.9) 

The differential operators in (2.8), (2.9) are called shift operators because of their param- 
eter shifting property. The differential operator in (2.8) followed by the one in (2.9) yields 
the second order differential operator of which the Jacobi polynomials are eigenfunctions. 
If we would have defined Jacobi polynomials only by their orthogonality property, not by 
their explicit expression, then we would have already been able to derive (2.8), (2.9) up 
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to constant factors just by the remarks that the operators D- in (2.8) and D^'^ in (2.9) 
satisfy 

f{D_f){x) g{x) (1 - x)- +1 (1 + xf +1 dx 

= - J* f(x) (D^g)(x) (1 - x)« (1 + xf dx (2.10) 

and that D- sends polynomials of degree n to polynomials of degree n — 1, while 
sends polynomials of degree n — 1 to polynomials of degree n. 

The same idea can be applied again and again for the more general classical orthogonal 
polynomials we will discuss in this section. 

It follows from (2.8), (2.9) and (2.10) that 



J (P^\x)f (1 - x) a (1 + x) p dx = const. J (1 - x) a+n (1 + xf +n 



dx, 



where the constant on the right hand side can be easily computed. What is left for computa- 
tion is a beta integral, which of course is elementary. However, we emphasize this reduction 
of computation of quadratic norms of classical orthogonal polynomials to computation of 
the integral of the weight function with shifted parameter, because this phenomenon will 
also return in the generalizations. Usually, the computation of the integral of the weight 
function is the only nontrivial part. On the other hand, if we have some deep evaluation 
of a definite integral with positive integrand, then it is worth to explore the polynomials 
being orthogonal with respect to the weight function given by the integrand. 

2.2. The Askey tableau. Similar to the classification discussed in §2.1, one can classify 
all systems of orthogonal polynomials {p n } which are eigenfunctions of a second order 
difference operator: 

a(x)p n (x + 1) + b(x)p n (x) + c(x)p n (x - 1) = X n p n (x). (2.11) 

Here the definition of orthogonal polynomials is relaxed somewhat. We include the pos- 
sibility that the degree n of the polynomials p n only takes the values n = 0, 1, . . . , N and 
that the orthogonality is with respect to a positive measure having support on a set of 
N + 1 points. Hahn [16] studied the (/-analogue of this classification (cf. §2.5) and he 
pointed out how the polynomials satisfying (2.11) come out as limit cases for q f 1 of his 
classification. 

The generic case for this classification is given by the Hahn polynomials 



Q n {x;a,p,N) := 3 F 2 



—n,n + a + f3+l,—x 
a+l,-N 



(-n)fc (n + a + g + l) fc {-x) k 
^ (n+l),(-.V ),,/,•! • 1 • 



k=0 
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Here we assume n = 0, 1, . . . , A and we assume the notational convention that 



TP 
r- 1 s 



-n, a,2, . . . , a r 



. \- (-n) fc (a 2 )fc...(ar)fc k „ a y / 91 ^ 
^ (&i)fc...(& a ) fc fc! 



remains well-defined when some of the bi are possibly non-positive integer but < —n. For 
the notation of Hahn polynomials and other families to be discussed in this subsection we 
keep to the notation of Askey & Wilson [8, Appendix] and Labelle's poster [24]. There 
one can also find further references. 

Hahn polynomials satisfy orthogonality relations 



N 1 

y^Qn{x)Q m {x)p{x) = 5 n , m — , n,m = 0,1..., A, (2.14) 

' — ' 7T. 



x=0 



where 



and 



PV ^ W (« + /3 + 2) JV V ; 

/A\ 2n + q + g + 1 (a + l) n (a + g + l) n 
7Fn ' W « + /3 + l (/3 + l) n (A + a + /3 + 2) n ' l ' J 

We get positive weights p(x) (or weights of fixed sign) if a, j3 e (—1, oo) U (— oo, —A). 

The other orthogonal polynomials coming out of this classification are limit cases of 
Hahn polynomials. We get the following families. 

Krawtchouk polynomials 

K n (x;p,N) := 2 F 1 (-n,-x;-N;p- 1 ), n = 0, 1, . . . , A, < p < 1, (2.17) 

with orthogonality measure having weights x i— > (^) p x (1 — p) N ~ x on {0, 1, ... , A}; 
Meixner polynomials 

M n (x;(3,c):= 2 F 1 (-n,-x;(3;l-c~ 1 ), < c < 1, /3 > 0, (2.18) 

with orthogonality measure having weights x i— > (,9) x c x /x! on Z + ; 
Charlier polynomials 

C n (x;a) := 2 Fo(-n,-x;-;-a~ 1 ), a > 0, 

with orthogonality measure having weights a; i— > a x /x! on Z+. 

Note that Krawtchouk polynomials are Meixner polynomials (2.18) with (3 := q~ N . 
Krawtchouk and Meixner polynomials are limits of Hahn polynomials, while Charlier poly- 
nomials are limits of both Krawtchouk and Meixner polynomials. In a certain sense, the 
very classical orthogonal polynomials are also contained in the class discussed here, since 
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Jacobi, Laguerre and Hermite polynomials are limits of Halm, Meixner and Charlier poly- 
nomials, respectively. For instance, 

P^ f3 \l-2x) = ( a+ p n lim Q n (Nx;a,p,N). 

See Table 1 (or rather a part of it) for a pictorial representation of these families and their 
limit transitions. 

The Krawtchouk, Meixner and Charlier polynomials are self-dual, i.e., they satisfy 

p n (x) =p x (n), x,n e Z + or x,n G {0, 1, . . . , N}. 

Thus the orthogonality relations and dual orthogonality relations of these polynomials 
essentially coincide and the second order difference equation (2.11) becomes the three 
term recurrence relation after interchange of x and n. Then it can be arranged that the 
eigenvalue X n in (2.11) becomes n. By the self-duality, the L 2 completeness of the systems 
in case N = oo is also clear. 

For Hahn polynomials we have no self-duality, so the dual orthogonal system will be 
different. Observe that, in (2.12), we can write 

k-i 

(-n)fc (n + a + (3 + l) fc = f[ (-n(n + a + (3 + 1) + j(j + a + (3 + 1)), 

j=0 

which is a polynomial of degree k in n(n + a + j3 + 1) . Now define 

R n (x(x + a + f3 + 1); a, /3, N) := Q x (n; o>, [3 , N) , n,x — 0,1,..., iV. 
Then R n (n = 0, 1, . . . , N) extends to a polynomial of degree n: 

—n, —x, x + a + (3 + 1 



:1 

a + l,-N 



R n (x(x + a + (3 + 1); a, f3, N) = 3 F 2 



(a+l) k (-N) k k\ 

These are called the dual Hahn polynomials. The dual orthogonality relations implied by 
(2.14) are the orthogonality relations for the dual Hahn polynomials: 

N 1 
R m {x{x + a + (3 + 1)) R n {x{x + a + P + l)) 

Here tt x and p(n) are as in (2.15) and (2.16). Thus the dual Hahn polynomials are also 
orthogonal polynomials. The three term recurrence relation for Hahn polynomials trans- 
lates as a second order difference equation which is a slight generalization of (2.11). It has 
the form 

a(x)p n (X(x+l)) + b(x)p n (X(x)) + c(x)p n (X(x-l)) = X nPn (X(x)), (2.19) 
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where A(x) := x(x + a + (3 + 1) is a quadratic function of x. The Krawtchouk and Meixner 
polynomials are also limit cases of the dual Hahn polynomials. 

It is a natural question to ask for other orthogonal polynomials being eigenfunctions 
of a second order difference equation of the form (2.19). A four-parameter family with this 
property are the Racah polynomials, essentially known for a long time to the physicists as 
Racah coefficients, which occur in connection with threefold tensor products of irreducible 
representations of the group 577(2). However, it was not recognized before the late seventies 
(cf. Wilson [32]) that orthogonal polynomials are hidden in the Racah coefficients. Racah 
polynomials are defined by 



R n (x(x + 7 + 5 + 1); a, /3, 7, S) := 4 F 3 



—n, n + a + j3 + 1, —x, x + 7 + 5 + 1 
a + 1,/3 + 5 + 1,7 + 1 



1 



(2.20) 



where a + lor/3 + 5 + lor7 + l = —N for some N e Z+, and where n = 0, 1, . . . , N. 
Similarly as for the dual Hahn polynomials it can be seen that the Racah polynomial R n 
is indeed a polynomial of degree n in A(x) := x(x + 7 + 5 + 1). The Racah polynomials 
are orthogonal with respect to weights w(x) on the points X(x), x = 0, 1 . . . , N, given by 



w(x) 



(7 + 5 + l) x ((7 + 5 + 3)/2) x (q + l) x ((3 + 5 + l) x (7 + 1), 
x\ ((7 + 5 + l)/2) x ( 7 + 5 -« + 1)^(7-/3 + 1), (<? + !)* 



It is evident from (2.20) that dual Racah polynomials are again Racah polynomials with 
a, /5 interchanged with 7, 5. Hahn and dual Hahn polynomials can be obtained as limit 
cases of Racah polynomials. 

Each orthogonality relation for the Racah polynomials is an explicit evaluation of a 
finite sum. It is possible to interpret this finite sum as a sum of residues coming from a 
contour integral in the complex plane. This contour integral can also be considered and 
evaluated for values of N not necessarily in Z+. For suitable values of the parameters the 
contour integral can then be deformed to an integral over the imaginary axis and it gives 
rise to the orthogonality relations for the Wilson polynomials (cf. Wilson [32]) defined by 



W n (x 2 ;a,b, c, d) 
:= (a + b) n (a + c) n (a + d) n 4-F3 



—n, n + a + b + c + d — 1, a + ix,a — ix 
a + b,a + c,a + d 



;i 



(2.21) 



Apparently, the right hand side defines a polynomial of degree n in x 2 . If a,b,c,d have 
positive real parts and complex parameters appear in conjugate pairs then the functions 
x 1— > W n (x 2 ) (n G Z + ) of (2.21) are orthogonal with respect to the measure w{x)dx on 
[0, cx)), where 



to(x) 



r(a + ix) r(6 + ix) r(c + ix) r(d + ix) 



r(2zx) 



The normalization in (2.21) is such that the Wilson polynomials are symmetric in their 
four parameters a, 6, c, d. 

The Wilson polynomials satisfy an eigenfunction equation of the form 



a(x) W n ((x + i) 2 ) + fe(x) W n (x 2 ) + c(x) W n ((x - i) 2 ) = A n W n (x 2 ). 
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So we have the new phenomenon that the difference operator at the left hand side shifts 
into the complex plane, out of the real interval on which the Wilson polynomials are 
orthogonal. This difference operator can be factorized as a product of two shift operators 
of similar type. They have properties and applications analogous to the shift operators 
for Jacobi polynomials discussed at the end of §2.1. They also reduce the evaluation of 
the quadratic norms to the evaluation of the integral of the weight function, but this last 
problem is now much less trivial than in the Jacobi case. 

Now, we can descend from the Wilson polynomials by limit transitions, just as we did 
from the Racah polynomials. On the 3F2 level we thus get continuous analogues of the 
Hahn and dual Hahn polynomials as follows. 

Continuous dual Hahn polynomials: 



S n (x 2 ; a, 6, c) := (a + b) n (a + c) n 3F2 



—n, a + ix, a — ix 
a + 6, a + c 



: 1 



where a, 6, c have positive real parts; if one of these parameters is not real then one of the 
other parameters is its complex conjugate. The functions x 1— > S n (x ) are orthogonal with 
respect to the measure w(x)dx on [0, 00), where 



w(x) := 



r(a + ix) F(b + ix) F(c + ix) 



T(2ix) 



Continuous Hahn polynomials: 



( 7 — T\ - n (a + a) n (a + b) n 
p n {x;a,b,a,b) :=i ; 3^2 



ni 



—n, n + a + a + b + b— l,a + ix 
a + a, a + b 



;i 



where a, b have positive real part. The polynomials p n are orthogonal on R with respect to 
the measure \T(a + ix) T(b + ix)\ 2 dx. In Askey & Wilson [8, Appendix] only the symmetric 
case (a, b > or a = b) of these polynomials occurs. The general case was discovered by 
Atakishiyev & Suslov [9]. 

Jacobi polynomials are limit cases of continuous Hahn polynomials and also directly of 
Wilson polynomials (with one pair of complex conjugate parameters). There is one further 
class of orthogonal polynomials on the 2-F1 level occurring as limit cases of orthogonal 
polynomials on the 3F2 level: 

Meixner-Pollaczek polynomials: 



Pi?\x;<l>) := 



(2a) 



2 F 1 (-n,a + ix; 2a; 1 - e 2i(p ), a > 0, < 4> < n. 



(Here we have chosen the normalization as in Labelle [24], which is the same as in Pol- 
laczek's original 1950 paper.) They are orthogonal on R with respect to the weight function 
x I— > e( 2< ^ _7r ) x |T(a + ix)| 2 . They can be considered as continuous analogues of the Meixner 
and Krawtchouk polynomials. They are limits of both continuous Hahn and continuous 
dual Hahn polynomials. Laguerre polynomials are limit cases of Meixner-Pollaczek poly- 
nomials. Note that the last three families are analytic continuations, both in x and in 
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the parameters, of dual Hahn polynomials, Halm polynomials and Meixner polynomials, 
respectively. 

All families of orthogonal polynomials discussed until now, together with the limit 
transitions between them, form the Askey tableau (or scheme or chart) of hypergeometric 
orthogonal polynomials. See Askey & Wilson [8, Appendix], Labelle [24] or Table 1. See 
also [19] for group theoretic interpretations. 

2.3. Big g-Jacobi polynomials. These polynomials were hinted at by Hahn [16] and 
explicitly introduced by Andrews & Askey [4]. Here we will show how their basic properties 
can be derived from a suitable pair of shift operators. We keep the convention of §1 that 
< q < 1. 

First we introduce (/-integration by parts. This will involve backward and forward 
q-derivatives: 

(D-f)(x) ■= f{x) - f{qx) (D+f)(x) ■= f(q ~ lx) - f(x) 

w q sm- {l _ q)x , Kv q j)w- {1 _ q)x ■ 

Here D~ coincides with D q introduced in (1.9). 

Proposition 2.1 If / and g are continuous on [— d, c] (c, d > 0) then 



ijqX. 



f {D~f){x) g(x) d q x = /(c) giq-'c) - f(-d) gi-q^d) - I f(x) (D+g)(x) d q 
J-d J-d 

Proof 

oo 

(D-f)(x) g(x) d q (x) = £(/(cg*) ~ f(cq k+1 )) g(cq k ) 

fc=0 

= limj /(c) g(q-i C ) - f(cq N+1 ) g(cq N ) + £ f(cq k ) { 9 (cq k ) - g{cq k ~ 1 )) } 
^ fc=0 ' 

oo 

= /(c) giq-'c) - /(0) g(0) + £ f(cq k ) {g(cq k ) - (/(eg*" 1 )) 

k=0 

= f(c)g(q- 1 c)-f(0)g(0)- [ f(x)(D+g)(x)d q x. 

Jo 

Now apply (1.13). □ 
Let 

/ . . x {qx/c.-qx/d-.q)^ 

w{x;a,b,c,d;q) := . (2.22) 

(qax/c, -qbx/d; q)^ 

Note that w(x; a, 6, c, d; q) > on [—d, c] if c, d > and 

< a < - & - — <b <- or [a = ca & b = da for some a E C\R.] (2.23) 
dq q cq q J 
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From now on we assume that these inequalities hold. If / is continuous on [—d, c] then 

lim / f(x)w(x;q a ,qP,c,d;q)d q x= [ f (x) (I - c' 1 x) a (I + d' 1 xf dx , (2.24) 
«ti J-d J-d 

Thus the measure w(x;a,b,c,d;q)d q x can be considered as a (/-analogue of the Jacobi 
polynomial orthogonality measure shifted to an arbitrary finite interval containing 0. 
Since 

w(q~ 1 c; a, 6, c, d; q) = = w(—q~ 1 d; a, 6, c, <i; q), 
we get from Proposition 2.1 that for any two polynomials f,g the following holds. 

c 

D qf)( x ) g( x ) w ( x ; Q a , c > ^ q) d q x 



q 

d 



f(x) [D+ (g(.) w(.; qa, qb, c, d; q))] (x) d q x. (2.25) 

d 



Define 

_ [ D g g Q > g & > c ? ^ g) /(•))] ( x ) 



(D+> a > b f)(x) 



w(x; a, b, c, d; q) 
--(1 - q)' 1 x _1 (l - x/c) (1 + x/d) /(g~V) 

- (1 - q)' 1 x~ l (1 - qax/c) (1 + qbx/d) f(x). (2.26) 



In the following D q f(x) orD"(/(x)) will mean (D~ f)(x). Also Dp a > b f(x) or Dp a > b (f(x)) 
will mean (D+> a ' b f)(x). A simple computation yields the following. 



= "J— ^ > ( 2 - 27 ) 

q~ n ci~^ — (/& 

= — — r-z — (qax/c; q) n + terms of degree < n — 1 in x. (2.28) 

[i. q)a 

Define the big q- Jacobi polynomial P n (x; a, 6, c, <i; g) as the monic orthogonal polyno- 
mial of degree n in x with respect to the measure w(x; a, 6, c, d; (/) d q x on [— d, c]. Later we 
will introduce another normalization for these polynomials and we will then write them as 
P n instead of P n . It follows from (2.24) that 



limP n (x; q a , q^ \ c, d; q) = const. P^ ( 



d — c + 2x 
d + c 



a Jacobi polynomial shifted to the interval [— d, c]. Big (/-Jacobi polynomials with d = 0, 
c = 1 are called little q-Jacobi polynomials. 

Two simple consequences of the definition are: 

P n (-x; a, 6, c, d; q) = (-l) n P n (x; 6, a, d, c; g) (2.29) 

and 
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P n (Xx; a, 6, c, d; q) = X n P n (x; a, b, A _1 c, A _1 <i; g), A > 0. 

It follows from (2.25), (2.27) and (2.28) that D~ and Dp a > b act as shift operators on the 
big g-Jacobi polynomials: 



D q P n (x; a, 6, c, d; q) = P n -i{x; qa, qb, c, d; q), (2.30) 



1-g 

D+' a > b P n _! (x; qa, qb, c, d; q) = _ H P n (x; a, 6, c, d; g). (2.31) 



Composition of (2.30) and (2.31) yields a second order (/-difference equation for the big 
g-Jacobi polynomials: 



D+> a >" D~ P n (x; a, 6, c, d; ? ) = g(1 ^3^^ a, 6, c, d; 9 ). (2.32) 

The left hand side of (2.32) can be rewritten as 

A(x) ((D-) 2 P n )(q- 1 x) + B(x) (D-Pn^q^x) 

= a(x)P n (q- 1 x) + b{x)P n {x) + c(x)P n (qx) (2.33) 

for certain polynomials A, B and a, 6, a Here 

A{x) = q~ x (1 — qax/c) (1 + qbx/d), 

_ (1 - qb)c - (1 - ga)d - (1 - q 2 ab)x 
B{X) ~ (l-q)cd ' 1 ] 

Compare (2.32), (2.33) with (2.1), (2.11). Apparently we have here another extension of 
the concept of classical orthogonal polynomials. Two other properties point into the same 
direction. First, by (2.30) the polynomials D~P n+ i (n = 0, 1,2, . . .) form again a system 
of orthogonal polynomials. Second, iteration of (2.31) yields a Rodrigues type formula. 
It follows from (2.26) that 

\1 a J qad(± — q) \QO>J 

Combination with (2.31) yields the recurrence 



q 2 ab — q n+1 ~ ( c 

P n — ;a,b, c, d;q 



(1 — q)cd \qa 



(l-qa)(l + qad/c) ~ ( c 

Pn-i ^r;qa,qb,c,d;q . 



qad(l — q) \q 2 a 
By iteration we get an evaluation of the big g- Jacob i polynomial at a special point: 

P n (± i a,t,cd iq ) = (±)'' (^f : (^f^ . (2 .35) 
\qa J \qaj {q n+1 ab;q) n 
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Now we will normalize the big g-Jacobi polynomials such that they take the value 1 at 

c/(qa): 

P n (x; a, 6, c, d; q) 



P n (x;a,b,c,d;q) := 

P n (c/(qa);a,b,c,d; q) 

Note that the value at —d/(qb) now follows from (2.29) and (2.35): 



(2.36) 



( -d 



P n — r;a,6,c,rf;g = — — 



\qb 



It follows from (2.29) and (2.37) that 



ad\ n (qb;q) n (-qbc/d;q) 7 



be J (qa;q) n (-qad/c; q) n ' 



(2.37) 



P n (-x; a, b, c, d;q)= ( J q \ n - qbc / d ' q ^ n p n ( x - b, a, d, c; «?). 



6c / (qa;q) n (-qad/c;q) r 

The following lemma, proved by use of (1.10), gives the g-analogue of a Taylor series 
expansion. 



Lemma 2.2 If f(x) := J2 k = Q c k (Q ax / c i a ) k then 



(W'K?^)^'- 1 )' (f)\™^ 



(?;?)/ 



ff) 



Now put f(x) := P n (x;a,b,c,d;q) in Lemma 2.2 and substitute (2.30) (iterated) 
and (2.36). Then the c k can be found explicitly and we obtain a representation by a 
g-hypergeometric series: 



P n (x;a,b,c,d; q) = Y^ 



k=0 



= 3<P2 



(q n ;q) k (q n+1 ab;q) k (qax/c;q) k q k 
(qa;q) k (-qad/c; q) k (q;q) k 

q~ n , q n+1 ab, qax/c 
qa, —qad/c 



;q,q 



(2.38) 



Andrews & Askey [4, (3,28)] use the notation Pn a,f3 \x;c,d : q), which coincides, up to a 
constant factor, with our p n (x; q a , q@ , c, d; q). 

It follows by combination of (2.29) and (2.38) that 



P n (x;a, 6, c, d;q) 
P n (-d/(qb);a,b,c,d; q) 



3<P2 



q n ,q n+ ab,-qbx/d 
qb, —qbc/d 



(2.39) 



where the denominator of the left hand side is explicitly given by (2.37). 

Next we will compute the quadratic norms. By (2.25), (2.30) and (2.31) we get the 
recurrence 

/ c ~ q n ~ 1 (l — q n )cd f c ~ 

(P%w)(x;a,b,c,d;q)d q x = — - _ qn+1(xh J (P n -i w )( x l ? a » qb, c, d; q) d q x 



g n(n - 1)/2 (g;g)n( cd) 

(q" l ab:q)„ 



n re 



w(x;q n a,q n b,c,d;q) d q x, (2.40) 
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where the second equality follows by iteration. Now the (/-integral of w from —d to c can 
be rewritten as a sum of two 20i's of argument q of the form of the left hand side of (1.42). 
Evaluation by (1.42) yields 

f° ( i a \ a n \ {q,-d/c,-qc/d,q 2 ab;q) 00 

/ w{x;a,b,c,d;q)d q x = (1 -q)c- t— — — . (2.41) 

J-d \1 a i Qo, -qbc/d, -qad/c; q)^ 

Together with (2.40) this yields: 

fl d (PZw)(x;a, b, c, d;q) d q x = ^ (n _ 1} /2 (q, go, qb, -qbc/d, -qad/c; q) n 

f^ d w(x;a,b,c,d;q)d q x ' ' (q 2 ab;q) 2n {q n+1 ab;q) n 

Now we can compute the coefficients in the three term recurrence relation for the big 
g-Jacobi polynomials P n (x) := P n (x; a, 6, c, d; q): 

xP n (x)=P n+1 (x) + B n P n (x)+C n P n - 1 (x). (2.43) 

Then C n is the quotient of the right hand sides of (2.42) for degree n and for degree n — 1, 
respectively, so 

_q n - 1 (l-q n )(l-q n a)(l-q n b)(l-q n ab)(d+q n bc)(c + q n ad) 
n ~ (1 - q 2n ~ l ab) (1 - q 2n ab) 2 (1 - q 2n + l ab) ' ( ' ' 

Then we obtain B n by substution of x := c/(qa) in (2.43), in view of (2.35). 

2.4. Little g-Jacobi polynomials. These polynomials are the most straightforward 
g-analogues of the Jacobi polynomials. They were first observed by Hahn [16] and studied 
in more detail by Andrews & Askey [3]. Here we will study them as a special case of the 
big q- Jacobi polynomials. 

When we specialize the big g-Jacobi polynomials (2.36) to c = 1, d = and normalize 
them such that they take the value 1 at then we obtain the little q-Jacobi polynomials 

P n (x;b,a,l,0;q) 

p n (x; a, 6; q) := — — — r 

P„(0;6,a,l,0;g) 

=2Ml~ n ^ n+1 ab;qa;q,qx) (2.45) 



=(- ? &)-» 9 - n(n - 1)/2 fo^sfc 
(?a;?)n 



g n ,q n+1 ab,qbx 
g6,0 



(2.46) 



Here (2.45) follows by letting d -> in (2.39), while (2.46) follows from (2.38) and (2.37). 
With the equality of (2.45) and (2.46) we have reobtained the transformation formula 
(1.38) for terminating 2<Pi series. From (2.46) we obtain 

Pn {q-'b-\a^q) = (-qb)- q -»W {qh]q)n 



(qa;q) n ' 
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From (2.45) and (1.37) we obtain 



p n (l;a,b;q) := (-a) n q^ 1 ^ (<?M) ' 



(qa;q) n ' 

Little (/-Jacobi polynomials satisfy the orthogonality relations 

1 f\ u. „a „8. _n„ /,._« „fl. _wa (9*5?) 



S,(a +1,0 + 1) 7 



jf p n (t; ? ) p m (t; <A q) t a Jffi^ <V 



. (2.47) 



The orthogonality measure is the measure of the g-beta integral (1.27). For positivity and 
convergence we require that 0<a<l,6<l (after we have replaced q a by a and q^ by 6 in 
(2.47)). It is maybe not immediately seen that (2.47) is a limit case of (2.42). However, we 
can establish this by observing the weak convergence as d J, of the normalized measure 
const. xw(x; q@, q a , l,d;q)d q x on [—d, 1] to the normalized measure in (2.47) on [0,1]: 

lim SU^ +lt ^)^(t;q^q a ,l,d;q)d q t = ^ {q^- q) n (-q^d; q) n 



dio £ d w{t;qt 3 ,q a ,l,d;q)d q t d i° (q a+f3+2 ;q)n 

~ (<r*+ 2 ; q)n - B q (a + 1, f} + 1) J {Q *' Q)n (qP+H; qU 

Here the first equality follows from (2.41) and (2.22), while the last equality follows from 
(1.27). 

Little (/-Jacobi polynomials are the (/-analogues of the Jacobi polynomials shifted to 
the interval [0,1]: 

lrm Pn (x;q ,^. ?) = 

If we fix b < 1 (possibly or negative) and put a = q a then little (/-Jacobi polynomials 
tend for q f 1 to Laguerre polynomials: 

limpn [ \ \ ; q a , b; q ' 



7n' n \ i-b "* '"'V l«(o)" 

In particular, little (/-Jacobi polynomials p n (x; a, 0; q), called Wall polynomials, are q- 
analogues of Laguerre polynomials. 

Analogous to the quadratic tansformations for Jacobi polynomials (cf. [14, 10.9(4), 
(21) and (22)]) we can find quadratic transformations between little and big (/-Jacobi 
polynomials: 

( 2. —1 2. 2\ 

P 2n (x; a, a, 1, 1; q) = ^l^^q^) ' ^ 
d / it \ xp n (x 2 ;q,a 2 ;q 2 ) 

P^^a,!,!;^ q . (2.49) 
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The proof is also analogous: by use of the orthogonality properties and normalization of 
the polynomials. 

Just as Jacobi polynomials tend to Bessel functions by 



2 

ar 



P K/3)( 1 ) V '4iV 2 , 

►Vi(-; a + 1; -(Ax/2) 2 ) = (Ax/2)" Q r(a + 1) J a {\x), n N /N -> A as JV -> oo, 



little Jacobi polynomials tend to the Hahn-Exton (/-Bessel functions (1.60) 

,n+l 



lim p JV _ n (g iV x;a,6;g)= lim 20i(? _iV+n , a6^- n+1 ; ag; g, g JV+1 x) 

iV— >oo N^oo 



= i0i (0; ag; g, g n+ x), 
cf. Koornwinder & Swarttouw [23, Prop. A.l]. 

2.5. Hahn's classification. Hahn [16] classified all families of orthogonal polynomials 
p n (n = 0, 1, . . . , N or n = 0, 1, . . .) which are eigenfunctions of a second order g-difference 
equation, i.e., 

A(x) (D 2 p n )(g" 1 x) + B{x) {Dqp^q-^x) = X n Pn(x) (2.50) 

or equivalently 

a(x)p n (q~ 1 x) + b(x)p n (x) + c(x)p n (qx) = \ n p n (x), 

where A, B and a, 6, c are fixed polynomials. Necessarily, A is of degree < 2 and B is 
of degree < 1. The eigenvalues A n will be completely determined by A and B. One 
distinguishes cases depending on the degrees of A and B and the situation of the zeros 
of A. For each case one finds a family of g-hypergeo metric polynomials satisfying (2.50) 
which, moreover, satisfies an explicit three term recurrence relation 

xp n (x) =p n+1 (x) + B n p n (x) +C n p n - 1 (x). 

(Here we assumed the p n to be monic.) Then we will have orthogonal polynomials with 
respect to a positive orthogonality measure iff C n > for all n. A next problem is to 
find the explicit orthogonality measure. For a given family of g-hypergeometric polyno- 
mials depending on parameters the type of this measure may vary with the parameters. 
Finally the limit transitions between the various families of orthogonal polynomials can be 
examined. 

In essence this program has been worked out by Hahn [16], but his paper is somewhat 
sketchy in details. Unfortunately, there is no later publication, where the details have all 
been filled in. In Table 2 we give a q-Hahn tableau: a g-analogue of that part of the Askey 
tableau (Table 1) which is dominated by the Hahn polynomials. In the r <p s formulas in the 
Table we have omitted the last but one parameter denoting the base except when this is 
different from q. The arrows denote limit transitions. Below we will give a brief discussion 
of each case. We do not claim completeness. 
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Several new phenomena occur here: 

1. Within one class of g-hypergeo metric polynomials we may obtain, depending on the 
values of the parameters, either a (/-analogue of the Jacobi-Laguerre-Hermite class or of 
the (discrete) Hahn-Krawtchouk-Meixner-Charlier class. 

2. q- Analogues of Jacobi, Laguerre and Hermite polynomials occur in a "little" version 
(corresponding to Jacobi polynomials on [0, c], Laguerre polynomials on [0, oo) and Hermite 
polynomials which are even or odd functions) and a "big" version ((/-analogues of Jacobi, 
Laguerre and Hermite polynomials of arbitrarily shifted argument). 

3. There is a duality under the transformation q \— > q , denoted in Table 2 by dashed 
lines. We insist on the convention < q < 1, but we can rewrite (/ _1 -hypergeometric 
orthogonal polynomials as (/-hypergeometric orthogonal polynomials and thus sometimes 
obtain another family. 

4. Cases may occur where the orthogonality measure is not uniquely determined. 

We now briefly discuss each case occuring in Table 2. 

1) Big q-Jacobi polynomials P n (x; a, 6, c, d; q), defined by (2.36), form the generic case in 
this classification. A(x) and B(x) in (2.50) are given by (2.34) and A n is as in the right 
hand side of (2.32). From the explicit expression for C n in (2.44) we get the values of 
a, 6, c, d for which there is a positive orthogonality measure. For c, d > these are given 
by (2.23). 

The q-Hahn polynomials can be obtained as special big q- Jacobi polynomials with 
—qad/c = q~ N , n = 0,1, ... ,N (N G Z+). For convenience we may take c = qa. The 
g-Hahn polynomials are usually (cf. [15, (7.2.21)]) notated as 



Q n (x;a, b.N;q) := 3 4>2 



q n ,abq n+1 ,x 

- N ; 9 

aq,q JV 



(g n ;q)k(abq n+1 ;q)k(x;q)k k 
^ (aq;q)k(q~ N ;q)k(q;q)k 



Here the convention regarding lower parameters q N (N E Z+) is similar to the convention 
for (2.13). 

g-Hahn polynomials satisfy orthogonality relations 

Y.{Q n Qm){q- x ; a, b, N; q) ^ g) « f N ~ x (aq)- x = 0, n * m. 

^ (q;q)x(q;q)N-x 

The weights are positive in one of the three following cases: (i) b < q~ l and < a < (/ -1 ; 
(ii) a,b > q~ N ; (iii) a < and b > q~ N . For q f 1 the polynomials tend to ordinary Hahn 
polynomials (2.12): 

limQ n (q- x ; q a , q? , N; q) = Q n (x; a, (3, N). 

2a) Little q-Jacobi polynomials p n (x; a, b; q), given by (2.45), (2.46), were discussed in §2.4. 
When we put b := q~ N_1 in (2.46), we obtain (/-analogues of Krawtchouk polynomials 
(2.17): the q-Krawtchouk polynomials 



q 71 , — b ^q n , x 
K n (x;b, N;q) := 3 2 ' _ N ' ;q,q 

v, q 



limQ n (x;a,-(qba) 1 ,N;q) 

a— >0 
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with orthogonality relations 

JV 



JV / -N \ 

J2(KnK m )(q- x ; 6, N; q) 19 ^ h {-b) x = 0, n + m, n, m = 0, 1, . 



.,N, 



x=0 

and limit transition 

liming-*; 6, N; q) = K n (x; 6/(6 + 1), N). 

See [15, Exercise 7.8(i)] and the reference given there. Note that the g-Krawtchouk poly- 
nomials K n (q~ x ; 6, N; q) are not self-dual under the interchange of x and n. 

2b) Big q-Laguerre polynomials, these are big g-Jacobi polynomials with 6 = 0: 



P n (x;a,0, c,d;q) = 3 2 



q n ,0, qax/c 
qa, —qad/c 
1 



(2.51) 



(-q- n ca- 1 d~ 1 ;q) 
Here the second equality follows by (1.39). Note that 



2<Pl 



q n ,c/x 
qa 



;q,-qx/d 



limP n (a;; q a , 0, c, (1 — q) 1 ; q) = i-Fi(— n; a + 1; c — x) = const. L"(c — x). 
'/'i 

Another family of (/-analogues of Krawtchouk polynomials, called afEne q-Krawtchouk 
polynomials, can be obtained from (2.51) by putting —qad/c = q~ N : 



K^f(x;a,N;q) := 3 </> 2 



q n ,0,x 



aq,q- 



-N 



Q n (x;a,0, N; q) 



with orthogonality relations 

f2( K n ff K^ f )(q- x ;a,N;q) 

x=0 

and limit transition 



(aq;q) x (aq) 

(q;q)x (q;q)N- 



= 0, n^m, n, m = 0, 1, . . . , N, 



]imK^f(q- x ; a, N; q) = K n (x; l-a,N). 



an 



See [15, Exercise 7.11] and the references given there. 
2c) q-Meixner polynomials 



M n (x;a,c;q) := 2 0i 



q n 7 x — q n+1 



(/a c 



(2.52) 
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For certain values of the parameters these can be considered as (/-analogues of the Meixner 
polynomials (2.18), cf. [15, Exercise 7.12]. However, when we write these polynomials as 



M n (qax;a,a ;q) = 2 <Pi 



q ,qax n+1 _ 

;q,-q a 

qa 



lim P n (x; a, ad, l,d;q) 



and a is complex but not real, then these polynomials also become orthogonal (not doc- 
umented in the literature) and they can be considered as (/-analogues of Laguerre poly- 
nomials of shifted argument. Moreover, these polynomials can be obtained from the big 
g-Laguerre polynomials (2.51) by the transformation q i— > q~ x . Therefore, we call these 
polynomials also big (/ _1 -Laguerre polynomials. 

For a := q~ N_1 the polynomials (2.52) become yet another family of (/-analogues 
of the Krawtchouk polynomials, which we will call aflme q~ x -Krawtchouk polynomials, 
because they are obtained from afflne (/-Krawtchouk polynomials by changing q into q~ x . 
These polynomials, written as 

M n (x; q- N ~\ -ft" 1 ; q) := 2 0i ((/"", x; q~ N ; q, bq n+1 ) = lim Q n (x; a, 6, N; q), 

a^oo 

where n = 0,l,...,iV, have orthogonality relations 

N (hn- n\ ( 1 \N-X „x{x-l)/2 

£(M B M m )(<T*; q -»~\ -ft" 1 ; q) l^p^Z*) I = 0, n * m, 

and limit transition 

limM n (q- x ; q- N ~\ -ft" 1 ; q) = K n (x; b~\ N). 

See Koornwinder [20]. 

3a) These orthogonal polynomials are not documented in the literature. 

3b) Wall polynomials (cf. Chihara [12, §VI.ll]) are special little (/-Jacobi polynomials 



p n (x;a, 0;q) = 2 </>i(q n ,0;qa;q,qx) 

= J -^ T ^,Mq- n ,x- 1 -,--,q,x/a). 



The second equality is a limit case of (1.39). Wall polynomials are (/-analogues of Laguerre 
polynomials on [0, oo), so they might be called little g-Laguerre polynomials. 

3c) Moak's [28] q-Laguerre polynomials (notation as in [15, Exercise 7.43]) are given by 
L«(x; q) := (g J + V )n i0i(<T n ; q a+1 ; q, -xq n+ « +1 ) 

Wi q)n 

1 2 Mq~ n ,-x;Q,q,q n+a+1 ). 



(q;q) 
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The second equality is a limit case of (1.38). They can be obtained from the Wall poly- 
nomials by replacing q by q~ l . Their orthogonality measure is not unique. For instance, 
there is a continuous orthogonality measure 

f°° x a dx 

/ L ^(x;q)K{x;q)—— r — = <J m>n , m ^ n, 

JO \~ l 1 ~~ Q) x 'i ?Joo 

but also discrete orthogonality measures 

r°° x a d x 

/ L m( CX '^) L n( CX '^) ~/ T\ \ ^ = S m,n, 171 ^ U, C> 0. 

Jo (-c(l - q)x; q)oo 

Chihara [12] Ch. VI, §2 calls these polynomials generalized Stieltjes-Wigert polynomials 
and he uses another notation. For certain parameter values these polynomials may be 
considered as (/-analogues of the Charlier polynomials, see [15, Exercise 7.13]. 

3d) Al-Salam-Carlitz polynomials 

U { n a \x) = U^(x;q) := (-1)" <f (- 1 )/ 2 a n 2 ^(^, x" 1 ; 0; q, qx/a) 
(cf. Al-Salam & Carlitz [1]) satisfy the three term recurrence relation 

x U£\x) = Ui%(x) + (1 + a) q n U<?\x) - a q^ 1 (1 - q n ) U^x). 

Thus they are orthogonal polynomials for a < 0. They can be expressed in terms of big 
(/-Jacobi polynomials by 

E#>(x) = P n (x;0,0,l,-a;g), 

so they are orthogonal with respect to the measure (qx, qx/a; q)ood q x on [a, 1], cf. (2.22). 
By the above recurrence relation these polynomials are (/-analogues of Hermite polynomials 
of shifted argument, so they may be considered as "big" (/-Hermite polynomials. 

3e) Al-Salam-Carlitz polynomials 

V^\x)=V^{x;q) ^uPfoq- 1 ) = (— l) n q~ n ^ n ~ 1 ^ 2 a n 2</>o((/ _n , x; — ; q n a~ 1 ) 

(cf. [1]). For a > these form another family of (/-analogues of the Charlier polynomials. 
On the other hand, the polynomials x i— > Vn (— qax) can be considered as (/-analogues 
of Hermite polynomials of shifted argument. 

4a) Stieltjes-Wigert polynomials 

S n (x; q) := q~< 2n+1 V 2 ^ {q~ n ; 0; q, -q n+s / 2 x) 

(cf. Chihara [12, §VI.2]) do not have a unique orthogonality measure. It was already 
noted by Stieltjes that the corresponding moments are an example of a non-determinate 
(Stieltjes) moment problem. After suitable scaling these polynomials tend to Hermite 
polynomials as q | 1. 
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4b) Al-Salam-Carlitz polynomials lli^ with a := —1: 

Ut'Hx) = q «»-W 2 Mq~ n , x- 1 ; 0; q, -qx) 
= P n (x;0,0,l,l;g) 

(cf. [1]), (/-analogues of Hermite polynomials. By the last equality there is a quadratic 
transformation between these polynomials and certain Wall polynomials. This is a limit 
case of the quadratic transformations (2.48) and (2.49). 

4c) Al-Salam-Carlitz polynomials Vn^ with a := —1: 

r n vt l) (ix) = r n <r n(n - 1)/2 2 <M<r n , ^; - ; -q n ) 



= lim P n (x; ic, ic, qc, qc, q) 

c^oo 

= x n 2 Mq- n ,q- n+1 ;0;q 2 ,-x- 2 ) 



(cf. [1]), also (/-analogues of Hermite polynomials. By the last equality there is a quadratic 
transformation between these polynomials and certain of Moak's (/-Laguerre polynomials. 
This is a limit case of (2.48) and (2.49). 

2.6. The Askey- Wilson integral. We remarked earlier that, whenever some nontrivial 
evaluation of an integral can be given, the orthogonal polynomials having the integrand as 
weight function may be worthwhile to study. In this subsection we will give an evaluation 
of the integral which corresponds to the Askey- Wilson polynomials. After the original eval- 
uation in Askey & Wilson [8] several easier approaches were given, cf. Gasper & Rahman 
[15, Ch. 6] and the references given there. Our proof below borrowed ideas from Kalnins 
& Miller [18] and Miller [26] but is still different from theirs. Fix < q < 1. Let 

^ / \ ( z i z '■> q)oo ^2 

(az,az~ 1 ,bz,bz~ 1 ,cz,cz~ 1 ,dz,dz~ 1 ;q) (X 

We want to evaluate the integral 

I r> (hz, 

I a ,b,c,d ■■= — <f> w a ,b,c,d{z) — , \a\, \b\, \c\, \d\ < 1. (2.54) 

2m J\z\ = l Z 

Note that I a ,b,c,d is analytic in the four complex variables a, b, c, d when these are bounded 
in absolute value by 1. It is also symmetric in a, b, c, d. 



Lemma 2.3 



1 — abed 

aM = (1 - ab) (1 - ac) (1 - ad) qaM ' 



(2.55) 
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Proof The integral 

W q 1 / 2 a,q 1 / 2 b,q 1 / 2 c,q 1 / 2 d( Z ) dz 



z — z 1 z 



i 

equals on the one hand 

/ W q 1 / 2 a,q 1 / 2 b,q 1 / 2 c,q 1 / 2 d(q 1/2z ) dz_ 
q x l 2 z — q~ 1 / 2 z~ 1 z 

f (1 - az) (1 - bz) (1 - eg) (1 - dz) w a , b}C}d (z) dz_ 
~ J q 1 / 2 z(l-z 2 ) z 

and on the other hand 

j 1 / 2 a,q 1 / 2 b,q 1 / 2 c,q 1 / 2 d( ( l~ 1/2z ) dz_ 

q~ 1 / 2 z — q 1 / 2 z~ 1 z 

(1 - a/z)(l - b/z)(l - c/z)(l - d/z)w a ,b,c,d( z ) dz 



w q , 



-/ 



gV2 z -i (1 - z -2) z 
Subtraction yields 







Q" 1/2 a _1 ~ abcd)(l - az)(l - az' 1 ) + (1 - ab)(l - ac)(l - ad)} 



dz 

xiflaAc,<i( z ) — = 2niq~ 1/2 a~ 1 {-(l - abcd)I qa ^ Cjd + (1 - ab)(l - ac)(l - ad)I aACjd }. □ 
By iteration of (2.55) and use of the analyticy and symmetry of I a ,b,c,d we obtain 

Iabcd= (abed-q)^ (256) 

"' , ' c ' (a&, ac, ad, be, bd, cd; q)oo ' ' ' 

We might evaluate /o,o,o,o by the Jacobi triple product identity (1.50), but it is easier to 
observe that I lq i/2_ 1 _ q i/2 = 1. (Note that this case can be continuously reached from 
the domain of definition of I a ,b,c,d in (2-54).) Hence (2.56) yields that /o,o,o,o = 2/(g; q)oo- 
Thus 

I 2(abcd; q)^ 

a ' ' c ' (ab, ac, ad, be, bd, cd, q; q)^ 

By the symmetry of w aj b, c ,d( z ) under z i— > z~ x we can rewrite (2.54), (2.57) as 



Here iu is still defined by (2.53). The integral (2.58) is known as the Askey-Wilson integral. 
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2.7. Askey- Wilson polynomials. We now look for an orthogonal system correspond- 
ing to the weight function in the Askey- Wilson integral. First observe that 



(az,az 1 ;q) k (bz,bz 1 ;q) l , dz 

w a ,b,c,d{ z ) — = 



^q k a,q l b,c,d 



2m J (ac,ad;q) k (bc,bd;q)i 



z (ac, ad; q) k (be, bd; q)i 
(ab;q) k+t 



(abed; q) k+l 



L a,b,c,d- 



Thus 



1 1 



2m 



' a,b,c,d 



E 

fc=0 



(g n ,q n 1 abcd;q) k q k (az,az 1 ;q) k \ (bz,bz 1 ;q) l 
(ab,q;q) k 



dz 



(ac,ad;q) k J (bc,bd;q)i Wa ' h ' c ' d ^ z 



= E 



(q n ,q n 1 abcd;q) k q k (ab;q) k+ i 



fc=0 

(ab; q)i 



(ab,q;q) k (abcd;q) k+ i 

q~ n , q n ~ 1 abcd, q l ab 
ab, q l abcd 



392 



(abed; q)i 
(ab;q)i (q~ n+1 c~ 1 d~ 1 , q~ l ; q) 



(abed; q)i (ab, q 



-n-l+l 



/(abcd);q) n 



q,q 

= 0, Z = 0,l,...,ra-1. 



(2.59) 



Here we used the g-Saalschiitz formula (1.61). 

The above orthogonality suggests to define Askey-Wilson polynomials (Askey & Wil- 
son [8]) 



p n (cos 9; a, b, c, d \ q) 
a~ n (ab, ac, ad; q) n 



=r n (cos 6; a, b,c,d\ q) 



E 

k=0 



q n ,q n 1 abed, ae %B , ae 
ab, ac, ad 

(q~ n , qlabcd; q) k q k (ae iB , ae~ 10 ; q) k 



(2.60) 
(2.61) 



(ab,q; q) k 



(ac, ad; q) k 



Since 



fc-i 



(ae lB ,ae- lB ;q) k = \\ (1 - 2aq J cos 6 + a 2 q 2j ), 
i=o 

formula (2.61) defines a polynomial of degree n in cos 9. It follows from (2.59) that the func- 
tions 6 I— > p n (cos6>; a, b,c,d \ q) are orthogonal with respect to the measure w a: b, c ,d(^ e ) d8 
on [0,7r], so we are really dealing with orthogonal polynomials. Now 

p n (x; a, b,c,d \ q) = k n x n + terms of lower degree, with k n := 2 n (q n ~ 1 abcd; q) n , (2.62) 

so the coefficient of x n is symmetric in a, b, c, d. Since the weight function is also symmetric 
in a, b, c, d, the Askey-Wilson polynomial will be itself symmetric in a, b, c, d. 
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Proposition 2.4 Let \a\, \b\, |c|, \d\ < 1. Then 

i r 

— / p n (cos 6) p m (cos 9) w (cos 9) d9 = 5 m ,nh n , 
27r Jo 



where 



p n (cos 9) = p n (cos 9; a, 6, c, d I q), 
w(cost') - 



and 



/i n (1 — q n ~ 1 abcd) (q, ab, ac, ad, be, bd, cd; q) n 
ho (1 — q 2n ~ 1 abcd) (abed; q) n 

(abed; q) 



h 



(q, ab, ac, ad, be, bd, cd; q) c 



The orthogonality measure is positive if a, b, c, d are real, or if complex, appear in conjugate 
pairs, 

Proof Apply (2.59) to 



T~ — ^~ f (PnPm) {(z + z 1 )/2;a,b,c,d \ q) w aACjd (z) 

*a,b,c,d 27TZ J 



dz 

z 



where p n is expanded according to (2.61) and p m similarly, but with a and b interchanged. 

□ 

Note that we can evaluate the Askey- Wilson polynomial p n ((z + z~ 1 )/2) for z = a 
(and by symmetry also for z = b, c, d): 

Pn(( a + a_1 )/2; a, b,c,d \ q) = a~ n (ab, ac, ad; q) n . (2.63) 

When we write the three term recurrence relation for the Askey- Wilson polynomial as 

2xp n (x) = A n p n+1 (x) + B n p n (x) + C n p n - 1 (x), (2.64) 

then the coefficients A n and C n can be computed from 

2/c n 2/c n _i h n 



kn+i k n h n —i 

where k n and h n are given by (2.62) and Proposition 2.4, respectively. Then B n can next 
be computed from (2.63) by substituting x := (a + a~ 1 )/2 in (2.64). 
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2.8. Various results. Here we collect without proof some further results about Askey- 
Wilson polynomials and their special cases and limit cases. 
The q-ultraspherical polynomials 



C n (cos0;fi | q) := Jjf , ,g ^ N ic>! 



q- n ,q n p 2 ,P 1 / 2 e ie ,p 1 / 2 e- i9 

^ 2 (q;q)n 

const. Pn (cos9;p 1 / 2 ,p 1 / 2 q 1 / 2 ,-p 1 / 2 ,-p 1 / 2 q 1 / 2 \ q) 

( I o\k i(n-2k)9 

(q;q)n t( ) (q 1 - n p-\q;q)k {q/p> 



are special Askey- Wilson polynomials which were already known to Rogers (1894), however 
not as orthogonal polynomials. See Askey & Ismail [5]. 

By easy arguments using analytic continuation, contour deformation and taking of 
residues it can be seen that Askey- Wilson polynomials for more general values of the 
parameters than in Proposition 2.4 become orthogonal with respect to a measure which 
contains both a continuous and discrete part (Askey & Wilson [8]) 

Proposition 2.5 Assume a, 6, c, d are real, or if complex, appear in conjugate pairs, and 
that the pairwise products of a, 6, c, d are not > 1, then 

1 f n 

— / p n (cOS 9) Pm(cOS 9) W (COS 9) d9 + y^p n (Xk)Pm(Xk)Wk = S m ,n h n , 
27T J Q V 

where ^(cos^), iy(cos6>) and h n are as in Proposition 2.4, while the Xk are the points 
(eq k + e~ 1 q~ k )/2 with e any of the parameters a, 6, c or d whose absolute value is larger 
than one, the sum is over the k G Z + with \eq k \ > 1 and Wk is Wk(ci',b,c,d) as defined 
by [8, (2.10)] when Xk = (aq k + a~ 1 q~ k )/2. (Be aware that (1 — aq 2h )/(l — a) should be 
replaced by (1 - a 2 q 2k )/(l - a 2 ) in [8, (2.10)].) 

Both the big and the little g-Jacobi polynomials can be obtained as limit cases of 
Askey- Wilson polynomials. In order to formulate these limits, let r n be as in (2.60). Then 

limr n ( 2 / ( ^^ I 9 

= P n (x;a, 6, c, d; q) 



and 



V ( q 1 X l/2\2 l/2\-2 1/2 1/2 1, I A (qb',q)n / , n, 

hmr n -r^-;q 1 A a,q 1 A , -q 1 , -q 1 b q ) = — — — — — p n (x; b, a; q . 

See Koornwinder [22, §6]. As A becomes smaller in these two limits, the number of mass 
points in the orthogonality measure grows, while the support of the continuous part shrinks. 
In the limit we have only infinitely many mass points and no continuous mass left. 

An important special class of Askey- Wilson polynomials are the Al-Salam-Chihara 
polynomials p n (x; a, 6, 0,0 | g), cf. [6, Ch. 3]. Both these polynomials and the continuous 
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g-ultraspherical polynomials have the continuous q-Hermite polynomials p n (x; 0, 0, 0, | q) 
(cf. Askey & Ismail [5]) limit case. 

The Askey- Wilson polynomials are eigenfunctions of a kind of g-difference operator. 
Write P n (e 10 ) for the expression in (2.61) and put 

_ (l-az)(l-bz)(l-cz)(l-dz) 
K) ' (1- z 2 )(l-qz 2 ) 

Then 

A(z)P n (qz) - (A(z) + Aiz-^P^z) + Aiz-^P^q-'z) 

= -(1 - q- n )(l - q n - 1 abcd)P n {z). 

The operator on the left hand side can be factorized as a product of two shift operators. 
See Askey & Wilson [8, §5]. 

An analytic continuation of the Askey- Wilson polynomials gives q-Racah polynomials 

'q- n q n+1 a[3,q- x ,q x+1 ~f5 



R n (q- x + q x+ ^5) := 4 03 



<xq,poq,iq 



where one of aq, f35q or ^q is q N for some A G Z + and n = 0, 1, . . . , N (see Askey & 
Wilson [7]). These satisfy an orthogonality of the form 

N 

^2 Rn{v{x)) R m (p(x)) w(x) = 5m jn h n , m, n = 0, 1, . . . , N, 

where p(x) := q~ x +q x+1 "f5. They are eigenfunctions of a second order difference operator 
in x. 

There is nice characterization theorem of Leonard [25] for the g-Racah polynomials. 
Let N G Z_|_ or N = oo. Let the polynomials p n (n G Z+, n < N + 1) be orthogonal with 
respect to weights on distinct points p k (k G Z + , k < N + 1). Let the polynomials p* n be 
similarly orthogonal with respect to weights on distinct points p* k . Suppose that the two 
systems are dual in the sense that 

Pn(^k) =P* k (p* n ). 

Then the p n are g-Racah polynomials or one of their limit cases. 
Exercises to §2 

2.1 Show that (at least formally) the generating function 

Hn{x) 
n=0 n! 

for Hermite polynomials (cf. [14, 10.13(19)]) follows from the generating function 



oo 

e 2xt-t 2 _ ST^ %W 



(1 - t)- a ~ X e - te /( 1 -*) = J2 Ln(x) t U , \t\ < 1, 

n=0 

for Laguerre polynomials (cf. [14, 10.12(17)]) by the limit transition 

H n (x) = (-l) n 2 n / 2 n\ lim oT n l 2 L^((2a) 1 / 2 x + a) 

a^oo 

given in (2.7). 
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2.2 Show that (at least formally) the above generating function for Laguerre polynomials 
follows from the generating function for Jacobi polynomials 

oo 

^Pj?>®{x)t n = 2 a+/3 R-^l-t + R^il + t + R)-?, \t\ < 1, 

n=0 

where 

R:= (l-2xt + t 2 ) 1/2 , 
(cf. [14, 10.8 (29)]) by the limit transition 



L${x) = lim P^\l-2x/(3) 

p— >oo 



given in (2.5). 



2.3 Show that (at least formally) the above generating function for Hermite polynomials 
follows from the specialization a = (3 of the above generating function for Jacobi polyno- 
mials by the limit transition 

H n (x) = 2 n nl lim a~ n / 2 P^ a \a~ l/2 x) 

given in (2.6). 

2.4 Prove that the Charlier polynomials 

C n (x;a) := 2 Fo(-n, -x; -; -a -1 ), a > 0, 
satisfy the recurrence relation 

xC n (x; a) = -aC n+1 (x; a) + (n + a) C n (x; a) - nC n _i(i; a). 

2.5 Prove that 

lim (-(2a) 1/2 ) n C n ((2a) 1/2 x + a) = H n (x) 

a^oo 

by using the above recurrence relation for Charlier polynomials and the recurrence relation 

H n+1 (x) -2xH n (x) + 2nH n _ 1 (x) =0 

for Hermite polynomials. 

2.6 Prove the following generating function for Al-Salam-Chihara polynomials (notation 
as for Askey- Wilson polynomials in (2.60), (2.61)): 

(zc, zd; g)oo _ z m Pmjcos 6>; 0, 0,c,d\ q) M<1 

2.7 Use the above generating function in order to derive the following transformation for- 
mula from little g-Jacobi polynomials (notation and definition by (2.45)) to Askey- Wilson 
polynomials by means of a summation kernel involving Al-Salam-Chihara polynomials: 

a n p n (cos 9; a, 6, c, d \ q) {ac^ad^q)^ 
(ab,ac,ad;q) n (ae i9 , ae~ ie ; q) ^ 

^ , m _i , _i , s m p m (cos9;0,0,c,d\ q) 
= ;q ab,q cd;q)a r . 
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Table 1. The Askey tableau 



201 (q 



2a) little q- Jacobi 
& (/-Krawtchouk 

< n+1 ab 



1) big (/-Jacobi & (/-Hahn 



q"'~ r± ab, qax; ga, 0; q) 
n+1 ab; qb; qx) 



2b) big (/-Laguerre 
& afhne (/-Krawtchouk 

> 2 (g-",0,^;ga,^;g) 
201 (q n , f ;ga; -gx/d) 



2c) big g _1 -Laguerre 
& (/-Meixner 
& afhne q~ l -Krawtchouk 

2 



i(q n , x; ga; —q n+1 j c) 




2</>l 



3a) 

x , q n b; 0; gx) 



3b) Wall 

2<Pi(q~ n ,0;qb;qx) 
2 0o(«?" n ,x" 1 ; ;x/6) 



3c) Moak & (/-Charlier 

i0i(g _n ; go; — q n+1 ax) 
(g" n ,-x;0;g n+1 a) 



2</>H 



3d) Al-Salam- 
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;0;gx/a) 
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